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FOREWORD 


In  order  Chat  photovoltaic  solar  cells  may  be  used  more  effectively  to 
provide  electrical  power  for  spacecraft,  it  is  desirable  to  refine  the  engineering 
design  of  solar  cells  and  power  systems ,  based  on  a  better  understanding  of  the 
space  radiation  environment  and  of  radiation  damage  effects  on  semiconductor 
components .  The  objective  of  this  handbook  is  to  provide  a  summary  of  some  of 
the  useful  analytical  methods  and  test  data  which  can  be  applied  in  designing 
radiation  resistant  power  systems .  It  is  found  that  the  most  serious  obstacle 
which  prevents  accurate  prediction  of  solar  cell  degradation  for  earth  satellites 
is  our  inadequate  knowledge  of  the  fluxes  and  energy  spectra  of  electrons  and 
protons  in  the  magnetosphere .  When  the  space  environment  can  be  better  defined , 
it  will  be  possible  to  predict  more  accurately  the  degradation  of  present  day 
types  of  silicon  solar  cells ,  for  which  the  radiation  damage  characteristics 
have  been  quite  well  determined  by  laboratory  research .  As  improved  types  of 
solar  cells  and  materials  become  available ,  it  will  be  necessary  to  make  laboratory 
measurements  of  their  radiation  resistance  to  protons  and  electrons  in  order  to 
predict  their  performance  in  space .  Therefore ,  the  performance  data  presented 
herein  may  become  obsolete .  However ,  it  is  hoped  that  the  methods  for  data 
correlation  and  analysis  presented  here  will  be  an  aid  to  understanding  the  nature 
of  the  problems  involved ,  and  will  be  useful  in  the  conduct  of  engineering  analysis . 
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I  INTRODUCTION 


Bombardment  of  solar  cells  by  energetic  electrons  and  protons  In  space  produces 
radiation  damage  which  decreases  their  power  output.  Earth  satellites  which  operate 
within  the  range  of  altitudes  between  about  400  and  40, 000  miles  may  encounter  appreci¬ 
able  fluxes  of  high  energy  charged  particles  which  are  trapped  in  die  earth's  magnetic 
field.  In  the  design  of  solar  cell  power  systems,  consideration  should  be  given  to  the 
electron  and  proton  fluxes  and  their  energy  spectra  along  the  proposed  orbit,  die  type 
of  solar  cells,  the  material  and  thickness  of  shielding  material,  the  solar  cell  operating 
voltage,  and  the  expected  operating  temperature  range.  This  handbook  is  intended  to 
provide  a  review  of  research  data  and  analytical  methods  which  can  be  used  to  desigp 
radiation  resistant  silicon  solar  cell  power  systems  for  earth  satellites  which  have 
orbits  passing  through  the  magnetosphere .  A  discussion  of  the  effects  of  solar  flares 
on  solar  cells  has  not  been  included.  For  a  discussion  of  this  subject,  one  may  refer 
to  a  report  by  Madey  (Ref.  53)  who  calculated  the  power  degradation  expected  from  a 
solar  flare. 

Appendices  are  included  which  summarize  data  on  radiation  darkening  of  solar 
cell  cover  slide  materials  and  radiation  effects  to  transistors  and  diodes.  Radiation 
effects  to  electronic  components  such  as  resistors  and  capacitors  are  not  diacuased 
because  the  radiation  damage  expected  in  the  space  environment  for  such  non-semi¬ 
conductor  components  is  generally  negligible. 


II  THEORY  OF  THE  SOLAR  CELL 


A.  The  Fhotovoltaic  Effect 

An  n-type  semiconductor  crystal  contains  a  small  concentration  of  doping  atoms 
which  ionize  and  supply  electrons  for  conduction.  A  p-type  crystal  contains  doping 
atoms  which  take  electrons  from  the  lattice,  leaving  "holes"  which  can  conduct  current 
by  moving  from  one  atom  to  the  next  as  an  adjacent  electron  moves  in. 

When  a  p-n  junction  is  produced  in  a  crystal,  electrons  from  the  n-side  thermally 
diffuse  into  the  p-region  where  they  are  called  minority  carriers.  Similarly,  holes  from 
the  p-side  become  minority  carriers  in  the  n-region.  These  minority  carriers  rapidly 
become  neutralized  by  recombining  with  the  majority  carriers.  Since  each  region  was 
initially  electrically  neutral,  the  loss  of  electrons  from  the  n-region  and  the  loss  of 
holes  from  the  p-region  each  contributes  to  establishing  a  potential  difference  across 
the  junction  with  a  plus  voltage  on  the  n-side.  This  built  in  electric  field  can  be  used 
for  photovoltaic  power  generation. 

When  a  photon  of  light  is  absorbed  in  a  crystal  it  will  ionize  an  atom  and  re¬ 
lease  an  electron,  thereby  producing  an  electron -hole  pair.  After  a  short  burst  of 
illumination,  the  electrons  and  holes  which  have  been  injected,  both  diffuse  until  they 
find  an  opposite  number  ar.d  recombine.  The  length  of  time  for  a  minority  carrier  to 
recombine  (the  lifetime)  depends  on  the  density  of  recombination  centers,  which  are 
crystal  defects  which  provide  sites  where  minority  carriers  axe  captured  and  then 
recombined  readily  with  majority  carriers. 

The  distance  travelled  by  a  minority  carrier  before  recombining  is  called  the 
minority  carrier  diffusion  length,  L,  which  is  related  to  the  lifetime,  t  ,  by: 

L  =  v'dT  (1) 

2 

where  D  is  die  diffusion  constant.  In  silicon,  the  value  of  D  is  about  38  cm  /sec.  for 
2  o 

electrons  and  13  cm  /  sec.  for  holes  at  300  K.  Both  L  and  t  are  dependent  on  the 
impurity  content  and  on  the  crystal  perfection,  which  is  affected  by  radiation  damage. 
They  also  are  functions  of  the  temperature. 
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When  light  is  absorbed  in  the  region  near  the  p-n  junction  of  a  solar  ceil, 
which  is  usually  0.25  to  1  micron  from  the  surface,  some  of  the  electrons  and  holes 
will  diffuse  to  the  junction.  At  this  point,  the  charges  will  be  separated  by  the  built 
in  electric  field  at  the  junction,  thereby  providing  a  current  which  can  flow  through 
an  external  load.  This  method  for  direct  conversion  of  photon  energy  to  electrical 
power  is  known  as  the  photovoltaic  effect. 

The  current  output  of  a  solar  cell  under  red  light,  which  produces  electron-hole 
pairs  several  hundred  microns  below  the  junction,  depends  on  the  magnitude  of  the 
diffusion  length  for  minority  carriers  in  the  base  region.  Therefore  exposure  to 
radiation  which  produces  recombination  centers,  primarily  affects  the  response  of  a 
solar  cell  to  red  light. 

B.  Solar  Cell  Design 

A  typical  solar  cell  consists  of  a  rectangular  wafer  of  high  purity  semiconductor 
crystal  with  a  p-n  junction  formed  near  the  front  surface  and  electrical  connections 
applied  to  the  front  layer  and  to  the  base  region,  as  shown  in  Fig.  1.  Usually  a  con¬ 
ducting  grid  is  used  on  the  front  surface  to  reduce  the  internal  resistance  to  lateral 
current  flow  in  the  thin  surface  layer. 

Until  recently  most  solar  cells  in  the  United  States  were  made  with  a  p-type 
layer  on  n-type  silicon.  The  development  of  radiation  resistant  n-on-p  silicon  cells 
by  Mandelkorn  (Ref.  54)  has  lead  to  their  selection  for  many  satellites. 

The  improved  performance  of  n-on-p  cells  under  radiation  can  be  attributed 
partly  to  the  fact  that  electrons,  which  are  the  minority  carriers  in  the  p-type  base 
material  have  about  a  three  times  greater  diffusion  constant  than  holes  do  in  the  n-type 
base  of  p-on-n  cells.  A  contr.butir.g  factor  is  that  the  types  of  recombination  centers 
produced  by  radiation  in  p-type  silicon  are  less  effective  in  shortening  the  minority 
carrier  diffusion  length  than  those  produced  in  n-type  silicon. 

Optimization  of  solar  cell  design  for  radiation  resistant  performance  requires 
locating  the  p-n  junction  only  about  0.25  micron  below  the  front  surface  to  maximize 
the  collection  of  current  generated  by  photons  near  the  blue  end  of  the  visible  spectrum, 
which  are  absorbed  and  produce  electron-hole  pairs  very  near  the  surface.  Cells  de¬ 
signed  in  this  way  are  called  shallow -diffused  or  blue -shifted  cells.  They  derive  a 
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smaller  fraction  of  their  power  from  the  red  end  of  the  solar  spectrum  and  are  there¬ 
fore  less  sensitive  to  degradation  of  the  diffusion  length  in  the  base. 

C.  Solar  Cell  Performance 

The  diode  equation  for  a  solar  cell  is: 


!  -  1,-1 

L  0  J 


where: 

I  =  current  through  load 
Iq  =  saturation  current  of  the  p-n  junction 
I  =  light  generated  current 
V  =  voltage  across  the  load 
q  =  electronic  charge  (esu) 

K  =  Boltzmann’s  constant 
A  =  a  constant  normally  between  1  and  3 
T  =  temperature  (°  K) 

This  equation  predicts  a  current  voltage  characteristic  as  shown  in  Fig.  2 

Under  short  circuit  conditions,  the  influence  of  1q  is  eliminated  and  the  abort 
circuit  current  is  equal  to  1^.  The  light -gene rated  current  is  a  function  of  the  light 
intensity,  the  absorption  coefficient  for  photons,  the  geometry  of  the  cell  and  the 
diffusion  length  for  minority  carriers  in  the  n  and  p  regions. 

Kleinman  (Ref.  48)  has  analyzed  the  effect  of  diffusion  length  on  the  short  circuit 
current.  His  theoretical  solution  indicates  a  variation  of  the  short  circuit  approximately 
proportional  to  the  logarithm  of  the  diffusion  length  over  the  range  from  L  *  10  to  200 
microns  which  is  of  practical  importance. 

It  has  been  found  experimentally  by  Space  Technology  Laboratories  (Ref.  28) 
that  data  on  short  circuit  current  as  measured  under  tungsten  light  can  be  correlated 
with  the  logarithm  of  diffusion  length  for  typical  p/n  and  n/p  cells.  These  data  are 
shown  in  Fig.  3A  and  3B 

Experiments  by  Bell  Telephone  Laboratories  have  yielded  data  which  permitted 
computation  of  the  short  circuit  current  under  space  sunlight  as  a  function  of  diffusion 
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length,  which  are  shown  in  Fig,  4  for  comparison  with  the  STL  correlation  line.  It 
is  noted  that  there  is  some  departure  from  linearity  with  log  L  for  the  BTL  data  on 
blue-shifted  p/n  and  n/p  cells . 

A  very  significant  factor  to  be  noted  in  Fig.  4  is  that  the  short  circuit  current 
does  not  decrease  as  rapidly  with  diffusion  length  reduction  under  space  sunlight  as 
when  the  measurement  is  made  under  tungsten  light ,  which  is  excessively  rich  in  red 
light. 

Madey  (Ref.  53)  has  used  BTL  data  (Ref.  75)  to  correlate  maximum  power  under 
space  sunlight  with  the  logarithm  of  diffusion  length,  as  shown  in  Fig.  5. 

Based  on  the  empirical  evidence,  it  is  assumed  that  the  prediction  of  solar  cell 
performance  in  a  radiation  environment  can  be  reduced  to  a  calculation  of  the  degradation 
of  diffusion  length ,  coupled  with  reference  to  data  which  show  how  the  performance 
parameters  of  the  particular  type  of  cell  in  space  sunlight  will  vary  with  diffusion  length . 

The  use  of  minority  carrier  diffusion  length  (or  lifetime)  to  correlate  radiation 
damage  effects  is  recommended  because  it  is  an  accurately  measurable  quantity  (see 
next  section)  which  is  sensitive  to  damage  by  penetrating  radiation  and  is  independent 
of  other  solar  cell  parameters  like  surface  optical  reflectivity,  surface  recombination 
velocity,  and  junction  depth,  all  of  which  influence  measurements  of  current  and 
power . 

In  nearly  all  space  applications  a  glass ,  silica  or  sapphire  cover  slide  is  used 
over  the  cell  to  aid  in  radiative  heat  rejection  and  to  carry  a  spectrally  selective  filter 
as  well  as  to  provide  radiation  shielding  and  protection  in  ground  handling.  Since  the 
efficiency  of  a  solar  cell  varies  (logarithmically)  with  light  intensity,  it  is  desirable  to 
use  a  cover  slide  material  and  adhesive  (if  necessary)  which  do  not  darken  unde  "  space 
radiation ,  or  to  estimate  the  decrease  in  transmittance  expected  in  the  wavelength  range 
of  interest  (from  0 .35  to  1.1  microns) .  Pertinent  test  data  are  given  in  Appendix  A . 
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FIG.  4  SHORT  CIRCUIT  CURRENT  DENSITY  VS.  DIFFUSION  LENGTH  FOR  TWO  LIGHT  SOURCES 


Minority  Carrier  Diffusion  Length,  L  (microns) 


Maximum  Power  Output ,  (milliwatts) 

FIG  .  5  MAXIMUM  POWER  OUTPUT  AS  A  FUNCTION 
OF  DIFFUSION  LENGTH  FOR  N/P  SOLAR 
CELL  IN  SPACE  SUNLIGHT 
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Ill  RADIATION  DAMAGE  TO  SOLAR  CELLS 


This  section  includes  a  discussion  of  the  physics  of  radiation  damage,  the 
effect  of  radiation  on  minority  carrier  diffusion  length,  analysis  of  the  degradation 
of  diffusion  length  in  space  using  laboratory  data  on  damage  coefficients,  and  shows 
how  the  solar  cell  parameters  such  as  short  circuit  current,  maximum  power  and 
current  at  a  fixed  voltage  can  be  correlated  with  the  diffusion  length. 

A  .  Production  of  Defects  and  Recombination  Centers  by  Radiation 

When  electrons  at  energies  greater  than  145  Kev  and  protons  at  energies 
greater  than  98  ev  bombard  a  silicon  crystal,  they  can  displace  an  atom  from  the 
crystal  lattice,  producing  a  lattice  vacancy  and  a  recoil  atom  which  comes  to  rest 
as  an  interstitial  atom .  With  high  energy  electrons  or  protons ,  the  recoil  atom  may 
have  enough  energy  to  displace  other  atoms  before  coming  to  rest.  The  displacement 
of  a  single  atom  is  called  a  Frenkel  defect .  The  displacement  of  two  adjacent  atoms 
is  called  a  divacancy .  When  very  high  energy  protons  (  500  Mev)  interact  with  a 
nucleus  they  can  cause  spallation  (star  production)  which  results  in  several  nuclear 
fragments  which  stop  at  short  range  in  the  crystal,  displacing  many  atoms . 

The  crystal  defects  produced  may  not  be  thermally  stable  except  at  very  low 
temperatures  .  At  room  temperature  the  vacancies  as  well  as  crystal  impurity  atoms 
can  migrate  and  form  stable  defects ,  some  of  which  act  as  recombination  centers  for 
electrons  and  holes  .  The  recombination  process  occurs  as  the  defect  center  first 
captures  a  minority  carrier  and  subsequently  captures  a  majority  carrier,  thereby 
annihilating  an  electron-hole  pair  .  The  silicon  A-center  is  a  stable  defect  consisting  of 
a  substitutional  oxygen  atom  occupying  a  site  in  the  silicon  lattice.  The  A-center  con¬ 
centration  is  dominant  in  determining  the  diffusion  length  of  minority  carriers  (electrons) 
in  p-type  silicon  when  it  is  irradiated  by  low  energy  electrons .  P-type  silicon  is  the  base 
material  in  radiation  resistant  n/p  solar  cells  and  less  susceptible  to  the  formation 
of  recombination  centers  than  n-type  silicon . 

The  silicon  E -center  appears  to  be  a  vacancy  next  to  a  substitutional  phosphorus 
atom  and  is  produced  predominantly  in  floating  zone  n-type  silicon .  Research  is 
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continuing  to  understand  the  nature  of  other  types  of  defect  centers  in  semiconductors 
and  to  determine  whether  it  is  possible  to  intentionally  dope  the  material  in  such  a 
way  that  the  stable  defects  do  not  act  as  recombination  centers  which  decrease  the 
diffusion  length  and  lifetime  (Ref .  79) . 

The  effectiveness  of  recombination  centers  in  shortening  the  diffusion  length 
in  silicon  is  generally  not  directly  proportional  to  the  number  of  atoms  which  are 
displaced  by  radiation,  but  depends  on  the  types  of  defects  that  are  introduced  by 
protons  and  electrons  of  different  energies,  which  affect  the  types  and  concentrations 
of  the  recombination  centers  that  are  formed  by  interaction  with  impurities  (Wysocki, 
Ref.  88,  has  shown  that  a  simpler  situation  prevails  in  gallium  arsenide,  where  the 
damage  produced  by  radiation  appears  to  be  proportional  to  the  number  of  atoms 
displaced .) 

The  rather  wide  scatter  in  the  susceptibility  of  various  types  ot  p/n  cells  to 
radiation  damage  is  believed  to  be  caused  by  inadequate  control  of  impurities  in  the 
n-type  base  material  which  affect  the  formation  of  recombination  centers. 

One  effective  way  to  improve  the  collection  efficiency  in  spite  of  the  presence  of 
recombination  centers  is  to  use  a  "drift-field"  solar  cell  in  which  the  doping  is  graded 
to  spread  out  the  electrostatic  field  region  at  the  junction  .  By  this  technique  (Ref.  20), 
minority  carriers  from  a  larger  volume  ot  the  cell  are  electrostatically  swept  to  the 
junction  and  are  less  susceptible  to  capture  by  the  recombination  centers  that  are 
present . 

B .  Simplifying  Assumptions 

An  assumption  which  is  ot  great  value  in  simplifying  the  analysis  of  solar 
cells  is  that  all  the  performance  parameters  of  a  particular  design  of  solar  cell  are 
unique  (unctions  of  the  minority  carrier  dilfusion  length  in  the  base  region.  This  is 
equivalent  to  the  assumption  that  the  only  effect  of  penetrating  radiation  on  a  solar  cell 
is  to  reduce  the  minority  carrier  diffusion  length  in  the  base  region  . 

This  basic  assumption  is'  found  to  be  valid  when  rhe  change  in  diffusion  length 
is  produced  by  penetrating  particles  which  generate  a  fairly  uniform  concentration  of 
defects  throughout  the  front  100  to  200  micron  thick  layer  of  a  solar  cell.  In  general, 
protons  above  a  few  Mev  in  energy  and  electrons  above  a  tew  hundred  K.ev  in  energy 
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impinging  on  bare  solar  cells  approximately  meet  this  requirement.  However,  low 
energy  protons  which  only  penetrate  a  few  microns  or  less  will  primarily  damage  the 
front  layer  and  junction  region  and  the  assumption  above  will  be  invalid.  This  more 
complicated  situation  is  not  treated  further  here  because  it  is  not  believed  to  be  of 
great  importance  when  cover  shields  are  used.  However,  it  is  pertinent  in  understanding 
the  behavior  of  bare  solar  cells  when  subjected  to  low  energy  protons  in  space  . 

C .  Effect  of  Radiation  on  Diffusion  Length 

The  measurement  of  diffusion  length  in  a  solar  cell  has  in  many  cases  been  made 
by  injecting  carriers  at  a  constant  low  rate  with  penetrating  ionizing  radiation  (e  .g. 
with  1  Mev  electrons  according  to  the  method  by  Rosenzweig,  Ref.  70)  and  measuring 
the  short  circuit  current  produced,  which  is  proportional  to  the  diffusion  length. 
Rosenzweig  used  a  12  mil  aluminum  foil  in  the  electron  beam  to  produce  fairly  uniform 
ionization  in  the  solar  cell .  A  reproducibility  of  measurements  within  +  3%  was 
achieved  using  1  Mev  electrons  . 

Measurements  of  diffusion  length  have  been  made  before ,  during  and  after 
bombardment  of  solar  cells  with  energetic  electrons  and  protons  in  particle  accelerators 
such  as  Van  de  Graaf  machines  and  cyclotrons  in  order  to  measure  radiation  damage . 

Experimentally  it  is  found  that  the  minority  carrier  diffusion  length  (L)  measured 
at  a  standard  temperature  varies  with  integrated  flux  of  penetrating  charged  particles 
(§)  according  to: 


(3) 


"0 


The  damage  coefficient  K  is  a  function  of  the  type  and  energy  of  the  particles  and  of  the 

material  and  its  impurity  concentrations .  K  is  by  definition  equal  to  d(l/L^)/dj>  and  is 

2  2 
the  change  in  ( 1/ 1.  )  introduced  per  unit  integrated  flux  (particle/cm  ) .  It  may  be  called 

the  damage  per  particle  (if  L  is  in  cm  instead  of  microns ,  K  has  the  units  of  particles  ) , 

Fig.  6A,  6B,  &  7  show  empirical  correlations  of  the  diffusion  length  for  various 

types  of  p/n  and  n/p  solar  cells  as  a  function  of  the  integrated  flux  of  1  Mev  electrons 

(Ref .  28) .  It  is  seen  that  after  initial  degradation  to  a  point  where  1  <  <  1 


"0 


the  data  can  be  correlated  by  the  approximate  equation:  — -  *=  K^>  or  L  -  (K$) 


-1/2 


(4) 
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FIG.frA  DIFFUSION  LENGTH  DEGRADATION  OF  N/P  SOLAR  CELLS 
WITH  1  MEV  ELECTRONS 


2 

Integrated  Flux  of  1  MEV  Electrons,  (E/cm  ) 

FIG.  6B  DIFFUSION  LENGTH  DEGRADATION  OF  P/N  SOLAR  CELLS 
WITH  1  MEV  ELECTRONS 
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FIG.  7  DEGRADATION  OF  DIFFUSION  LENGTH  IN  i  OHM-CM  P-TYPF,  SIUCON 
AS  A  FUNCTION  OF  INTEGRATED  FLUX 
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In  Fig .  6A  it  is  seen  that  the  data  for  n/p  cells  from  three  manufacturers 
using  1  ohm-cm  p-type  silicon  for  the  base  material  all  cluster  quite  closely  together. 
The  improved  radiation  resistance  of  cells  with  higher  base  resistivity  (higher  purity) 
is  also  clear,  with  one  possible  exception. 

In  Fig .  6B  it  is  noted  that  p/n  cells  exhibit  a  large  scatter  in  the  data ,  presumably 
because  of  variations  of  the  impurity  content  in  the  base  material.  (However,  Weller 
(Ref.  S3)  has  shown  that  oxygen  content  does  not  control  the  recombination  of  minority 
carriers  in  p/n  ceils  when  irradiated  with  4 .8  Mev  protons). 

D .  Electron  Damage  Coefficients 

The  electron  damage  coefficients  for  p/n  and  n/p  cells  under  Mev  electrons 
have  been  reported  by  Rosenzweig,  Gummel ,  and  Smits  (Ref.  71),  and  the  effect  of 
electron  energy  has  been  presented  by  Rosenzweig  (Ref.  72).  Damage  coefficients  for 
1  Mev  electrons  have  also  been  calculated  from  test  data  by  Denney  (Ref.  28)  and  Statler 
(Ref.  77) .  Fig.  8  and  9  summarize  the  data  available.  It  is  seen  in  Fig.  8  that  there 
is  a  close  verification  of  the  test  data  at  1  Mev  for  n/p  cells  of  several  base  resistivity 
values  among  the  various  sources .  However ,  Fig .  9  shows  that  there  is  a  wide  range 
of  Kg.  values  for  p/n  cells .  As  shown  by  Rosenzweig  (Ref.  72)  the  damage  coefficient 
for  p/n  cells  can  be  correlated  with  the  particular  type  of  silicon  used  as  die  base 
material .  He  has  also  compared  the  relative  damage  rates  at  various  electron  energies 
for  three  types  of  cells ,  as  shown  in  Fig .  10 .  These  data  indicate  that  electrons  above 
1  Mev  in  energy  do  an  anomalously  large  amount  of  damage  to  the  p-type  silicon  in 
n/p  cells .  This  effect  has  also  be  found  by  Denney  at  STL.  Wysocki  (Ref.  88)  has 
suggested  that  the  high  energy  electrons  (and  protons)  produce  divacancy  defects  which 
lead  to  a  more  deleterious  recombination  center  in  p-type  silicon  than  in  n-type . 

E  .  Proton  Damage  Coefficients 

Downing  and  Denney  (Ref.  29)  have  shown  evidence  that  the  measured  diffusion 
length  and  lifetime  apparently  increase  with  the  minority  carrier  density  (injection 
level),  for  material  irradiated  by  high  energy  protons .  Van  Lint  et  al  (Ref.  80)  have 
observed  that  this  phenomenon  does  not  seem  to  be  as  predominant  for  1  Mev  electron- 
irradiated  material  as  it  is  for  material  irradiated  by  high  energy  protons  or  higher 
energy  electrons  (which  produce  more  complex  defects).  The  question  is  not  yet  settled 
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as  to  whether  the  effect  of  injection  level  on  measured  diffusion  length  for 

proton  irradiated  specimens  is  dependent  on  the  energy  of  the  protons . 

A  technique  for  making  measurements  using  a  chopped  electron  beam  in 
combination  with  a  steady  light  source  to  increase  the  injection  level  has  been  developed 
by  Downing  and  Denney  (Ref .  29) .  However ,  the  fact  that  the  light  source  does  not 
inject  minority  carriers  uniformly  through  the  solar  cell  makes  intexpretation  of  the 
results  rather  difficult .  Except  where  noted  below,  the  values  of  proton  damage 
coefficients  have  all  been  based  on  measurements  of  diffusion  length  using  low  injection 
levels . 


Measurements  of  the  effect  of  degradation  of  diffusion  length  on  the  short  circuit 
current  and  maximum  power  of  solar  cells  under  proton  bombardment  have  been  made 
by  several  experimenters .  In  order  to  correlate  the  data  on  a  common  basis ,  it  was 
necessary  in  some  cases  to  compute  the  damage  coefficient  from  data  on  the 
degradation  of  short  circuit  current  or  of  maximum  power  obtained  under  a  variety  of 
light  sources .  (See  Section  III-G-3  for  the  correlation  method  used). 

Tables  lll-l  and  I II -2  show  the  resulting  values  of  the  proton  damage  coefficients 
at  various  proton  energies  for  n/p  andp  /n  cells ,  respectively .  These  data  have  been 
plotted  in  Fig.  11  and  12.  Despite  considerable  scatter  in  the  data,  approximate 
correlation  curves  have  been  drawn  under  the  assumption  that  the  damage  coefficient 
varies  inversely  with  energy  up  to  a  certain  energy  and  then  remains  constant .  This 
correlation  is  chosen  for  its  mathematical  simplicity  in  damage  calculations  and  may 
only  be  accurate  within  a  factor  of  two  or  three  in  the  range  from  20  to  60  Mev  which  is 
usually  of  considerable  importance  in  determining  proton  damage  rates  in  space.  The 
points  marked  with  vertical  arrows  which  fall  considerably  above  the  correlation  line 
may  be  in  error  because  the  diffusion  lengths  were  measured  at  low  injection  levels . 

By  using  data  by  Denney  (Ref.  27)  the  points  at  80  Mev  and  9S.S  Mev  can  be  corrected 
to  correspond  to  an  injection  level  near  that  produced  by  space  sunlight,  in  which  case 
they  fall  close  to  the  correlation  line . 

By  comparing  Fig.  11  and  12,  it  is  seen  that  the  correlation  lines  yield  a  ratio 
of  damage  coefficients  for  1  ohm -cm  n-type  silicon  (p/n  cells)  to  that  for  1  ohm -cm 
p-type  silicon  (n/p  cells)  of  four  (4)  at  proton  energies  ip  to  60  Mev . 


TABLE  111-1 

Proton  Damage  Coefficients  for  10  -cm  p-Type  Silicon  (n/p  Cells) 


Proton  Energy  MEV 


Damage  Coefficient, 


•V 


Reference 


1 .8 

3.2 x  10‘6 

6 

4  .8 

1.7  x  10'6 

83 

8.3 

2.1 x  10‘6 

7* ** 

16.8 

8.3  x  10‘7 

70 

17.6 

1 .0  x  10‘6 

6 

10 .0 

9  .2  x  10“7 

7* 

20.5 

4  .0  x  10  7 

23*' 

45.0 

1 .0  x  10  ^ 

57 

80 .0 

6.0  x  10"7 

57 

95.5 

2.4  x  IQ’7 

23 

95.5 

5.4  x  10‘7 

23 

120 

4. Ox  10"7 

57 

130 

3.3  x  10"7 

70 

450 

1 .8  x  10'7 

23 

750 

2.8  x  10"? 

23 

*  Kp  was  calculated  from  the  authors'  values  of  L^,  L  and  Q. 

**  The  0c's  from  six  types  of  cells  were  averaged  and  converted  to  a  K  value . 


TABLE  III- 2 

Proton  Damage  Coefficients  for  10  -cm  n-Type  Silicon  (p/n  Cells) 

Proton  Energy,  ME V  Damage  Coefficient,  Kp  Reference 


1.8 

1 .0  x  10"5 

6 

4.8 

8.7 x  10'6 

83 

8.1 

1 .2  x  10  ^ 

7* 

16.8 

5 .1  x  10  ^ 

70 

17.6 

5.9  x  10"6 

6 

20.8 

7.9  x  10’7 

23 

45.0 

6.0  x  10"6 

57 

80.0 

4  .0  x  10"6 

57 

95.5 

3.0  x  lO*7 

23 

95 .5 

2.26  x  10'6 

23 

130 

2.0  x  10‘6 

70 

450 

3.1  x  10'7 

23 

450 

1  .2  x  10'6 

57 

740 

4.2  x  10'7 

23 

•K 


P 


was  calculated  from  the 


authors'  values  of 


'O' 


L  and  <!> . 


Proton  Knerjjy  (MFV) 

I  K'..  II  I’ROIDN  DAMAGF  COF.  P  FIOIF.NT  FOR  N  T  Sll.I(X)N 
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FIG.  12  PROTON  DAMAGE  COEFFICIENT  FOR  P/N  SILICON 
CELLS  AS  A  FUNCTION  OF  PROTON  ENERGY 
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The  damage  coefficients  at  proton  energies  of  1  Mev  and  lower  are  only  of 
interest  for  bare  solar  cells,  which  degrade  very  rapidly  when  exposed  to  geomagnetically- 
trapped  protons . 

The  proton  damage  coefficients  as  a  function  of  proton  energy  E  may  be  approxi¬ 
mated  by  the  following  relations: 

For  1  ohm-cm  p-type  silicon  (n/p  cells): 

Kp  =  1 .5  x  10  ^  E  *  for  1  <  E  <  60  Mev 
Kp  =  2.5  x  10* 7  for  E  >  60  Mev 
For  1  ohm-cm  n-type  silicon  (p/n  cells): 

K  =  6  x  Ilf 5  E*1  for  1<  E  <  150  Mev 
V  -7 

Kp  =  4  x  10  for  E  >  150  Mev 

F  .  Calculation  of  the  Diffusion  Length  From  Laboratory  Data 

If  more  than  one  type  or  energy  of  bombarding  particles  is  present,  it  is  a 

2 

good  approximation  to  assume  that  the  damage  (tire  differential  change  in  1/L  )  is 
additive  and  independent  of  the  rate  of  irradiation  .  It  is  also  usually  assumed  that  the 
damage  is  independent  of  the  temperature  and  the  operating  condition  of  the  solar  cell 
during  irradiation ,  although  it  is  known  that  radiation  effects  can  anneal  out  at  elevated 
temperature . 

Under  these  assumptions,  when  several  types  of  particles  are  present 

simultaneously  or  in  sequence: 


“a  ‘  \  :  Ml  +  M2  + 


r 


(5) 


If  the  radiation  in  space  h  is  n  continuous  energy  spectrum,  it  is  convenient  to 
define  a  damage  integral  iy  the  eqi.  n: 

D=  y  i  lCvc,t)i:i:.)ci!  -1-ij  -  —  y>  (6) 


where  a(E,t)  =  ommdiiect.or.il  iiux  p  r  u  ii;  energy  interval  in  space 

2 

(particles/cm  “  -  sec  -Mev)  at  time  t  ,  but  ir.od.fied  by  the  shield  over  the  solar  cell. 
K(E)  =  damage  coefficient  at  energy  E  . 
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The  factor  of  1/2  is  introduced  to  account  for  the  usual  assumption  that  there 
is  infinite  shielding  on  the  back  side  of  the  solar  cell . 

A  laboratory  irradiation  to  a  time  integrated  flux  with  particles  of  energy 
E  is  equivalent  to  a  space  irradiation  if 

K  ,0,  =  \  i  [p(E  ,t)  K(E)  dEdt  (7) 

11  2  fc  Jo 


where  K  ^  is  die  damage  coefficient  for  particles  of  energy  E  . 

If  E  is  taken  as  1  Mev  electrons ,  Equation  7  defines  die  1  Mev  unidirectional 
normal  incident  integrated  flux  ^  which  is  equivalent  to  a  certain  irradiation  for  a 
time  t  in  space . 

For  shielded  solar  cells  ,  the  spectral  flux  p(E  ,t)  and  the  corresponding  damage 
coefficient  K  should  be  based  on  the  flux  and  energy  spectrum  present  behind  the  shield 
at  the  surface  of  the  cell ,  if  the  damage  coefficient  K  has  been  measured  for  bare  cells  . 
Also  it  is  assumed  that  the  particles  penetrating  the  shield  have  a  range  longer  than  about 
200  microns  and  produce  uniform  damage  in  the  sensitive  base  region  of  the  solar  cell. 

The  equivalence  between  doses  (integrated  fluxes)  of  particles  at  any  arbitrary 
energy  and  the  dose  at  energy  E^  can  be  established  by  the  relation: 


or: 


K,  <{>.  =K  $ 
11  n  Yn 


(8) 

(9) 


The  ratio  (K^/K  is  the  relative  damage  rate  for  particles  of  energy  En  as 
compared  to  energy  E  . 

G  .  Correlation  of  the  Performance  of  Solar  Cells  as  a  Function  of  Diffusion  Length 
1  .  Spectral  Response 

Fig.  13  shows  the  change  in  spectral  response  of  a  typical  n/p  1  ohm -cm  cell 
under  l  Mev  electron  bombardment.  It  is  seen  that  the  major  degradation  occurs  in  the 
response  near  the  red  end  of  the  spectrum  ,  where  photons  are  absorbed  more  deeply  in 
the  cell  and  minority  earners  must  diffuse  a  greater  distance  through  the  damaged 
material  to  reach  die  junction  .  The  value  of  the  minority  carrier  diffusion  lengtfi 
corresponding  to  each  value  of  integrated  flux  is  shown  on  the  figure  . 
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2 .  Current  -  Voltage  Characteristic 

Fig.  14  is  a  typical  plot  showing  the  change  in  the  current-voltage  (i-V)  character¬ 
istics  of  a  solar  cell  after  various  degrees  of  irradiation .  It  is  seen  that  the  short 
circuit  current  degrades  more  rapidly  than  the  open  circuit  voltage 
3  .  Short  Circuit  Current 

It  was  shown  in  Fig.  3A,  3B,  and  4  that  the  short  circuit  current  and  maximum 
power  of  solar  cells  can  be  correlated  with  the  diffusion  length .  However ,  the  particular 
correlation  obtained  depends  on  the  spectral  distribution  and  intensity  of  the  light  source 
used ,  and  the  presence  or  absence  of  an  anti-reflection  coating  on  the  cell  (which 
typically  increases  the  power  output  by  at  least  25%,  Ref.  75).  Also,  the  initial  value 
of  short  circuit  current  or  power  obtained  from  a  particular  cell  will  be  a  function  of 
the  initial  diffusion  length  of  the  particular  piece  of  silicon  from  which  it  was  made, 
as  well  as  of  differences  in  surface  recombination  properties  and  other  properties  due 
to  manufacturing  tolerances  . 

If  one  assumes  reasonably  close  control  of  the  manufacturing  process ,  then 
it  should  be  possible  to  select  an  average  value  for  the  initial  diffusion  length  and 
corresponding  average  values  of  the  initial  short  circuit  current  and  maximum  power 
as  measured  under  a  particular  light  source  which  are  representative  of  the  particular 
type  of  production  cells  .  It  is  then  desirable  to  predict  what  the  cell  performance  will 
initially  be  in  space  sunlight  and  how  the  performance  parameters  will  degrade  as  the 
diffusion  length  is  reduced  by  radiation . 

In  order  to  obtain  a  generally  applicable  procedure  for  predicting  performance , 
the  following  method  was  used: 

If  we  assume  that  short  circuit  current  (I)  for  a  particular  type  of  solar  cell 
varies  directly  with  the  logarithm  of  the  diffusion  length  (L),  then  one  can  write: 

1  =  C  In  (^  )  (6) 

L 

where  C  and  L  are  constants  which  depend  only  on  the  light  source. 

Then  the  short  circuit  current  ratio  is: 


Current  Density  (ma/cm 


which  can  also  be  written  as: 


4-  =  n-[ln<-T>]  <8> 

*0  1  L  1  L0  J 

This  equation  indicates  that  for  a  given  light  source,  which  determines  L', 
and  with  a  given  initial  diffusion  length  ,  the  term  In  (L^/L*  )  is  a  constant  and  the 
short  circuit  current  should  degrade  linearly  with  the  term  In  (L/L^). 

Fig.  15  shows  the  variation  of  short  circuit  current  ratio  with  the  logarithm 
of  the  ratio  of  final  to  initial  diffusion  length  for  three  different  light  sources .  The 
curves  shown  are  based  on  measurements  for  both  p/n  and  a/p  cells  by  Rosenaweig, 
Gummel  and  Smits  (Ref  .71),  Denney ,  Downing  and  Van  Atta  (Ref.  28)  and  Weller 
(Ref.  83). 

Data  from  Weller  (Ref.  83)  for  measurements  with  water-filtered  tungsten 
light  show  that  Equation  8  is  approximately  satisfied .  Data  computed  for  space  sunlight 
(Ref.  71)  appear  to  depart  somewhat  from  linearity  as  may  be  seen  in  Fig.  15  • 

However,  the  use  of  the  dimensionless  parameter  (L/L^)  to  correlate  test 
data  is  still  useful  as  a  technique  for  predicting  short  circuit  current  ratio .  Using 
the  curves  in  Fig.  15  it  is  only  necessary  to  know  to  predict  the  degradation  of 
short  circuit  current  ratio  with  L .  It  is  believed  that  the  curves  in  Fig.  15will  permit 
estimating  the  short  circuit  current  ratio  for  blue-shifted  silicon  solar  cells  correspond¬ 
ing  to  any  value  of  L  from  10  to  200  microns  with  reasonable  accuracy ,  independent 
of  whether  the  cell  is  n/p  or  p/n ,  with  base  resistivity  values  from  1  to  10  ohm-cm , 
and  with  mitial  diffusion  lengths  from  100  to  200  microns . 

The  problem  of  predicting  the  absolute  value  of  initial  short  circuit  current 
under  space  sunlight  from  ground  measurements  with  other  light  sources  is  not  con¬ 
sidered  here  .  (See  Ref .  71) 

By  cross -plotting  from  Fig.  15  ,  the  value  of  the  diffusion  length  required  to 
yield  a  25%  reduction  in  short  circuit  current  under  various  light  sources  is  plotted 
as  a  function  of  the  initial  diffusion  length  in  Fig.  16. 

2 

From  the  value  of  L  determined  from  Fig.  16.  the  change  in  (1/L  )  required  to 
produce  25%  degradation  in  short  circuit  current  under  various  light  sources  can  be 


60, 
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computed.  This  value  is  calculated  to  be  4.1x10  cm  based  on  data  for  blue  nfr 
cells  with  Lq  =  120 u  in  space  sunlight. 

4  .  Maximum  Power  Output 

A  typical  plot  of  the  variation  in  the  maximum  power  ratio  as  a  function  of  the 
short  circuit  current  ratio  is  shown  in  Fig.  n .  It  is  seen  that  the  maximum  power 
decreases  only  slightly  more  rapidly  than  the  short  circuit  current. 

Fig.  18  shows  the  power  output  of  n/p  solar  cells  as  a  function  of  cell  voltage 
at  various  degrees  of  degradation  produced  by  1  Mev  electrons  .  These  data  for  space 
sunlight  have  been  correlated  in  dimensionless  form  as  maximum  power  ratio  vs . 
diffusion  length  ratio  in  Fig.  19  .  along  with  the  comparable  data  for  filtered  and  un¬ 
filtered  2800°K  tungsten  light .  Data  computed  for  space  sunlight  (Ref.  71)  for  blue 
p/n  cells  with  on  initial  diffusion  length  of  ISO  u  are  included . 

5 .  Power  at  a  Fixed  Voltage 

The  rate  of  degradation  in  solar  cell  current  when  operated  at  a  fixed  cell  voltage 
(as  when  connected  to  charge  a  battery)  depends  on  the  value  of  the  operating  voltage  . 

If  the  operating  voltage  is  chosen  too  high ,  the  cell  current  will  degrade  very  rapidly . 

On  the  other  hand,  if  the  cell  voltage  is  too  low,  then  the  maximum  power  available 
may  not  be  obtained .  It  is  seen  from  Fig .  18  that  if  the  cell  is  to  operate  until  its 
diffusion  length  is  degraded  to  about  7  .6  microns ,  the  cell  voltage  should  be  selected 
at  approximately  0.40  volts  in  order  to  obtain  maximum  available  power  at  the  end  of 
life  .  However,  consideration  of  temperature  coefficients  may  force  selection  of  a 
lower  voltage . 

6 .  Temperature  Coefficients  of  Solar  Cells  as  a  Function  of  Radiation  Damage 

Martin  .  Teener  and  Ralph  (Ref.  55)  have  investigated  the  temperature  coefficient 

of  silicon  solar  cells  using  a  Spectrosun  solar  simulator .  They  found  that  the  temperature 
coefficient  of  short  circuit  current  for  n/p  10  ohm-cm  cells  increased  from  an  initial 
value  of  about  60  a  amp/°C  to  102 .5  uamp/°C  after  irradiation  with  a  flux  of  10^  1  Mev 
electrons/cm2  .  The  value  for  p/n  1  ohm-cm  cells  increased  from  62 .5  to  100  u  A/°C  . 
However .  the  temperature  coefficient  of  open  circuit  voltage  stayed  constant  at  about 
-2  .24  mv/°C  for  the  n/p  cells  and  -2 .05  mv/°C  for  the  p/n  cells  .  Ross  (Ref .  73) 
measured  a  temperature  coefficient  for  the  maximum  power  voltage  of  about  -2 .08  mV/°C 
for  both  n/p  and  p/n  cells  . 
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It  is  important  to  recognize  that  the  temperature  coefficient  of  short  circuit 
current  would  be  about  a  factor  of  ten  too  high  if  measured  with  a  tungsten  light 
source ,  which  overemphasizes  the  red  response  as  compared  to  sunlight  and  therefore 
overstresses  the  effect  of  variation  of  minority  carrier  lifetime  with  temperature . 

Unless  the  negative  temperature  coefficient  of  open  circuit  voltage  is  taken  into 
account,  the  operating  voltage  point  may  be  selected  too  close  to  the  open  circuit 
voltage  and  very  severe  degradation  of  performance  will  occur  as  radiation  damage 
proceeds . 

Experimental  results  on  the  combined  effects  of  temperature  and  radiation 
damage  on  the  1-V  characteristics  of  solar  cells  as  reported  by  Martin ,  Teener  and 
Ralph  (Ref.  55)  are  shown  in  Fig.  20,21,  22,  and  23.  By  reference  to  these  figures, 
the  operating  voltage  point  may  be  selected  so  that  maximum  power  is  obtained  after 
degradation  to  the  required  value  of  equivalent  1  Mev  electron  flux  (or  corresponding 
diffusion  length) . 


FIG.  20  CURRENT  VS.  VOLTAGE  CURVES  FOR  N/P  10  OHM-CM 

SILICON  SOLAR  CELLS  AT  VARIOUS  CELL  TEMPERATURES 
BEFORE  AND  AFTER  I  MEV  ELECTRON  IRRADIATION 
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FIG  .  21  CURRENT  VS  .  VOLTAGE  CURVES  FOR  P/N  1  OHM-CM 

SILICON  SOLAR  CELLS  AT  VARIOUS  CELL  TEMPERATURES 
BEFORE  AND  AFTER  1  MEV  ELECTRON  IRRADIATION 


-  41  - 


20  40  60  80  100 

CELL  TEMPERATURE  (  DEG.  C  ) 


FIG  .  22  ELECTRICAL  CHARACTERISTICS  AS  FUNCTIONS  OF 
TEMPERATURE  AT  INCREASING  LEVELS  OF  1  MEV 
ELECTRON  FLUX  FOR  N/P  SOLAR  CELLS 
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FIG.  23  ELECTRICAL  CHARACTERISTICS  AS  FUNCTIONS  OF 
TEMPERATURE  AT  INCREASING  LEVELS  OF  1  MEV 
ELECTRON  FLUX  FOR  P/N  SOLAR  CELLS 
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IV  THE  SPACE  RADIATION  ENVIRONMENT 

(See  also  recent  Reference  145  in  the  July  1963  issue  of  the  Bell  System  Technical  Journal) 


A.  Introduction 

The  most  damaging  space  radiations  for  solar  cells  are  the  geomagnetic&lly 
trapped  electrons  and  protons  in  die  magnetosphere  and  the  solar  flare  protons  which 
occur  sporadically  at  times  near  maximum  solar  activity  and  which  are  observed 
primarily  outside  the  magnetosphere .  Other  radiations  which  are  generally  negligible 
with  respect  to  radiation  damage  on  semiconductors  include  primary  galactic  cosmic 
rays  (which  are  primarily  protons  and  heavier  ions),  neutrons  produced  by  cosmic  ray 
interaction  with  the  atmosphere ,  solar  x-rays  and  gamma  rays . 

For  earth  satellites  within  the  magnetosphere,  it  is  important  to  know  the  flux  and 
energy  spectra  of  electrons  and  protons  encountered  in  order  to  calculate  radiation 
effects .  An  excellent  survey  article  as  of  October  1962  was  prepared  by  O'Brien 
(Ref .  104) .  Unfortunately,  recent  experimental  data  on  trapped  electrons  and  protons 
has  not  yet  been  properly  correlated  to  provide  a  good  knowledge  of  the  flux  and  energy 
spectrum  of  trapped  particles  at  all  locations .  Furthermore  the  "Starfish"  nuclear 
explosion  of  July  9 ,  1962  injected  large  numbers  of  artificial  electrons  which  are  still 
decaying  so  that  the  electron  fluxes  are  still  changing  with  time .  Therefore ,  the  brief 
summary  below  must  be  considered  tentative  and  subject  to  change  alter  further 
correlation  of  experimental  measurements . 

B.  Me  Ilwain's  Coordinates 

In  order  to  discuss  the  distribution  of  trapped  particles  it  is  most  convenient  to 
use  the  coordinate  system  developed  by  Me  Ilwain  (Ref.  99) .  The  two  coordinates  are 
B,  the  scalar  magnitude  cf  the  magnetic  field,  and  L,  a  parameter  that  is  very  nearly 
constant  along  each  magnetic  field  line  and  which  is  nearly  equal  to  the  radial  distance 
to  the  field  line  at  the  geomagnetic  equator.  Fig.  24  shows  an  example  of  L  lines  plotted 
with  respect  to  magnetic  dipole  coordinates  R  and  X ,  taken  from  (Ref.  99)  where  R 
is  the  distance  from  the  center  of  the  magnetic  dipole  and  X  is  the  geomagnetic  latitude . 
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FIG .  24  THE  B-L  MAGNETIC  COORDINATE  SYSTEM 


The  trapped  particles  may  be  visualized  as  moving  primarily  in  helical  paths 

around  each  magnetic  field  line  (  »  constant  L)  and  being  reflected  between  mirror 

points  north  and  south  of  the  magnetic  equator  at  points  where  the  field  B  has  increased 

2 

over  its  value  at  the  equator  by  the  ratio  (1/sin  c^)  where  orQ  IS  pitch  angle  of  particle  at 
magnetic  equator .  Each  particle  bounces  between  the  north  and  south  mirror  points 
with  a  period  of  the  order  of  seconds  or  less .  In  addition ,  electrons  drift  eastward  and 
protons  westward  with  a  period  of  the  order  of  a  minute  required  to  circle  the  earth, 
although  this  longitudinal  motion  can  be  neglected  in  understanding  the  steady  state 
flux  distributions  in  the  belt . 

The  fluxes  and  energy  spectra  of  particles  can  be  correlated  analytically  quite 
well  with  respect  to  the  B-  L  coordinate  system ,  The  pitch  angle  distribution  of  the 
particles  can  also  be  correlated  with  the  variation  in  the  magnitude  of  the  field  B  along 
each  line  of  constant  L .  The  flux  is  higher  perpendicular  to  the  magnetic  field  Mm  than 
parallel  to  the  field  line ,  particularly  at  the  mirroring  points .  This  effect  of  pitch  angle 
distribution  is  significant  for  a  satellite  (like  ANNA- IB)  which  is  oriented  with  respect 
to  the  magnetic  field .  However ,  generally  the  flux  may  be  defined  in  terms  of  an  average 
omnidirectional  flux  at  a  point  in  B-L  space . 


An  alternate  plot  of  the  relationship  between  the  B-L  coordinates  and  the  geomagnetic 
coordinates  R  and  X  is  shown  in  Fig.  25.  In  converting  from  B-L  coordinates  to  earth - 
centered  coordinates  it  must  be  recognized  that  the  earth's  magnetic  field,  if  approxi¬ 
mated  as  a  simple  dipole ,  is  centered  at  a  point  located  411 .4  Km  from  the  center  of 
the  earth  on  a  line  toward  the  point  at  150 .8  degrees  east  longitude  and  15.6  degrees 
north  latitude  (Ref .  3  ) .  The  axis  of  the  dipole  is  tilted  at  an  angle  of  11 .7  degrees 
from  the  earth’s  axis  toward  the  geographic  longitude  of  69  degrees  west . 

Even  if  fluxes  in  the  B-L  coordinates  were  known  precisely,  in  order  to  predict 
fluxes  at  an  arbitrary  point  in  space  with  an  accuracy  better  than  a  factor  of  10  it  is 
necessary  to  include  the  anomalous  variations  in  the  magnetic  field  of  the  earth  from 
the  dipole  approximation .  An  important  example  is  the  South  Atlantic  Anomaly,  where 
magnetic  field  lines  dip  to  low  altitudes  and  appreciable  fluxes  of  trapped  protons  and 
electrons  are  encountered  at  altitudes  of  only  1000  Km . 

C.  Electron  Fluxes  and  Energy  Spectra 

The  flux  of  natural  electrons  in  the  inner  region  of  the  magnetosphere  have  been 
augmented  by  artificial  electrons  injected  by  high  altitude  nuclear  explosions .  Table  IV- 1 
lists  the  explosions  which  have  produced  measurable  effects .  Of  these ,  only  die 
"Starfish"  explosion  on  July  9 ,  1962  has  produced  a  long-lasting  belt  of  electrons  which 
is  of  importance  in  affecting  solar  cells . 

Hess  (Ref.  96)  has  prepared  maps  of  die  artificial  electron  belt  in  B-L  and  R-X 
coordinates  (Fig. 26  and  21)  which  are  applicable  for  a  time  of  one  week  after  the 
'Starfish"  explosion  on  the  assumption  that  the  electrons  all  have  a  fission  spectrum 
as  shown  in  Fig. 28.  This  map  by  Hess  showed  appreciable  fluxes  of  electrons  at 
values  of  L  ig>  to  4  (altitudes  up  to  3  earth  radii  at  the  magnetic  equator) .  However , 
data  by  Van  Allen  (Ref.  Ill)  indicated  that  the  fission  electrons  were  confined  to  lower 
values  of  L .  It  now  appears  that  the  fission  spectrum  is  approximately  applicable  only 
for  L  less  than  2  and  that  the  electron  spectrum  is  softer  (a  steeper  variation  of  flux 
with  energy)  than  a  fission  spectrum  for  L  greater  than  2 ,  where  many  of  the  electrons 
observed  were  naturally  present  prior  to  "Starfish" .  An  effort  is  now  being  devoted  by 
Hess  at  the  Goddard  Space  Flight  Center  to  correlate  the  electron  spectra  as  a  function 
of  L  (Ref.  40) .  Until  this  correlation  becomes  available  it  seems  appropriate  to  assume 
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Plotted  according  to  the  transformation: 


R  =  L  cos  X 


FIG.  25  THE  MAPPING  OF  THE  POLAR  COORDINATES  R  and  X  ONTO 

THE  B  -  L  PLANE 
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FIG.  26  THE  B  -  L  MAP  OF  ELECTRON  FLUXES  APPROXIMATELY  ONE  WEEK 

AFTER  STARFISH 


that  Hess'  map  as  shown  in  Fig .  27  applies  approximately  for  values  of  L  less  than  2 
and  that  the  fission  spectrum  as  shown  in  Fig.  28  applies  in  this  region .  However,  it 
is  clear  that  the  electron  fluxes  have  continued  to  decay  since  one  week  after  Starfish 
and  that  the  flux  values  appropriate  at  a  time  of  5- 1/2  months  after  Starfish  (Dec .  1962) 
may  be  approximately  a  factor  of  two  lower  than  shown  in  Fig .  27  .  Until  further 
correlation  work  is  completed ,  it  does  not  appear  useful  to  try  to  calculate  electron 
fluxes  and  spectra  for  L  values  from  2  to  4 . 

The  average  electron  fluxes  in  the  outer  zone  (i.e. ,  for  values  of  L  greater 
than  4)  are  shown  in  Fig  .29  (Ref .  93) .  It  should  be  noted  that  the  electron  flux  in  this 
region  fluctuates  radically  during  each  day  Ref .  103)  and  also  as  a  result  of  geomagnetic 
storms  (Ref.  108)  and  that  the  applicable  energy  spectra  in  this  region  are  not  well 
established.  It  is  seen  in  Fig. 29  that  the  average  shape  of  the  constant  flux  contours  In 
the  magnetic  equatorial  plane  are  distorted  by  interaction  with  plasma  from  the  sun 
(the  solar  wind) . 

D.  Proton  Fluxes  and  Energy  Spectra 

Mcllwain  and  Pizzella  (Ref.  100)  have  shown  that  the  energy  spectrum  of 
protons  throughout  the  magnetosphere  may  be  approximated  by  the  equation: 

pp(E)dE  =  (Constant)  e‘E/E°dE  (9) 

with 

Eq  =  306  L  Mev  ( 10) 

This  equation  gives  a  reasonable  fit  to  the  data  for  values  of  L  up  to  8  and  is  in  not 
bad  agreement  with  Naugle  and  Kniffen’s  spectra  (Ref .  102)  obtained  with  an  emulsion 
experiment  at  L  values  from  1 .47  to  1 .79  . 

The  value  of  the  constant  which  defines  the  absolute  flux  at  any  point  may  be 
approximated  by  using  the  proton  flux  map  of  Fig.30(Ref.  101),  which  shows  contours 
of  the  integral  flux  above  31  Mev.  Improved  proton  maps  are  being  prepared  by  Walter 
L.  Brown  of  Bell  Telephone  Laboratories  and  by  the  Goddard  Space  Flight  Center  in 
cooperation  with  the  experimenters  . 

Additional  pertinent  references  on  the  space  environment  are  included  on  page  R-7. 
It  is  expected  that  improved  correlations  of  data  on  the  fluxes  and  energy  spectra  of 
trapped  radiation  will  become  available  in  the  near  future,  in  which  case  more  accurate 
calculation  of  space  radiation  effects  will  become  feasible . 
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FIG  29.  Summary  of  observed  omnidirectional  intensities  of  electrons  (E,J  40  kev) 
obtained  from  approximately  twenty  complete  orbits  of  Explorer  14  during 
October-December  1962  and  from  Explorer  12  during  August-December  1961. 
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FIG.  30  PROTON  FLUX  MAP 


TABLE  IV-1 

Data  On  Nuclear  Explosions  Which  Injected  Particles  Into  The  Magnetosphere 


Country 

Identification 

Size 

Location 

Date 

Remarks 

U.S. 

Teak  &  Orange 

1  MT 

Less  than  50  mi .  above 
Johnston  Island 

Aug  1  and 
12.  1958 

Effects  disappeared 
in  a  few  days  due 
to  particle 
absorption  in  the 
South  Atlantic 
Anomaly . 

U.  S. 

Argus  (3  shots) 

1-2 

KT 

300  mi .  above  South 
Atlantic 

1958 

Artificial  belts 
lasted  several 
weeks. 

U.S. 

Starfish 

1.4 

MT 

250  mi .  above 

Johnston  Island 

7/9/62 

Produced  long- 
lasting  belt  of 
artificial  electrons . 

USSR 

3  shots 

? 

High  altitude  and  high 
latitude 

10/22/62 

10/28/62 

11/1/62 

Artificial  electron 
belts  in  slot 
between  inner 
and  outer  zones 
decayed  rapidly . 

V  CORRELATION  OF  SATELLITE  TEST  DATA  ON 
SOLAR  CELL  PERFORMANCE 


This  section  summarizes  the  results  of  telemetered  performance  data  for 
solar  cells  with  various  thicknesses  of  transparent  shielding  on  satellites  which  have 
been  subjected  to  radiation  in  the  natural  and  artificial  radiation  belts. 

Table  V-l  shows  the  orbital  parameters  and  date  of  launch  of  fourteen  satellites 
for  which  telemetered  solar  cell  data  were  available.  In  most  cases,  the  telemetered 
data  were  available  in  the  form  of  short  circuit  current  with  normal  solar  incidence 
angle  as  a  function  of  time  in  orbit.  Although  the  temperature  of  the  solar  cells  was 
telemetered  in  some  cases,  no  effort  has  been  made  to  correct  all  data  to  a  standard 
temperature . 

An  example  of  the  variation  of  short  circuit  current  ratio  (I/I^)  with  time  for 
p/n  and  n/p  silicon  solar  cells  with  various  shield  thicknesses  on  three  satellites  is 
shown  in  Fig.  31.  The  flight  test  data  were  extrapolated  where  necessary  to  estimate 
the  time  in  orbit  required  for  the  short  circuit  current  to  decrease  by  25%.  These 
extrapolations  were  made  in  most  cases  by  assuming  a  20%  reduction  in  short  circuit 
current  per  decade  of  time  after  at  least  10%  degradation  had  occurred,  except  where 
test  data  showeda  different  trend.  This  extrapolation  procedure  is  inaccurate  when 
the  radiation  exposure  is  not  constant  with  time,  which  is  the  case  for  many  satellites 
where  apogee  precession  affects  the  flux  encountered  per  orbit.  Table  V-2  shows  the 
data  on  the  extrapolated  life  to  75%  of  initial  short  circuit  current  for  solar  cell3  on 
various  satellites.  The  shield  thicknesses  have  been  converted  to  grams  per  square 
centimeter. 

The  estimated  time  in  days  for  the  n/p  solar  cells  to  degrade  by  25%  have  been 
plotted  against  the  shield  thickness  in  grams  per  square  centimeter  in  Fig.  32.  A 
similar  plot  of  the  data  for  p/n  solar  cells  is  shown  in  Fig.  33. 

An  inspection  of  the  data  showed  that  for  the  Relay- 1  and  Anna  1-B  satellites, 
which  carried  both  n/p  and  p/n  cells,  there  was  a  factor  in  the  neighborhood  of  four 
times  longer  life  for  n/p  cells  than  for  p/n  cells  under  shield  thicknesses  greater  than 


TABLE  V-l 


ORBITAL  PARAMETERS  OF  SATELLITES 


Period 


No. 

Satellite 

Inclination 

Launch  Date 

I. 

Explorer  XI 

490.8  KM 
264  N.MI. 

1799  KM 

970  N.MI. 

28.8° 

Apr.  27,  1961 

2. 

Midas  III 

3450  KM 
1850  N.MI. 

3510  KM 

1890  N.MI. 

91.2° 

161.5 

July  12,  1961 

3. 

Explorer  XII 

304  KM 

163  N.MI. 

77, 000  KM 
41,500  N.MI. 

33° 

Aug.  16.  1961 

4. 

Midas  -  IV 

3530  KM 
1898  N.MI. 

3760  KM 

2025  N.MI. 

95.9° 

166 

Oct.  21,  1961 

5. 

TRAAC 

960  KM 

516  N.MI. 

1106  KM 

596  N.MI. 

32° 

105.6 

Nov.  15,  1961 

6. 

TRANSIT  -  IV  B 

960  KM 

516  N.MI. 

1106  KM 

596  N.MI. 

32° 

105.6 

Nov.  15.  1961 

7. 

Ariel  *  I 

390  KM 

210  N.MI. 

1210  KM 

650  N.MI. 

54° 

100.8 

Apr,  26,  1962 

8. 

Tfclstar  -  I 

952  KM 

513  N.MI. 

5660  KM 

3040  N.MI. 

45° 

157.7 

July  10.  1962 

9. 

Tetrahedral  112  N.MI. 

Research  Satellite  -  I 

362  N.MI. 

82° 

Sept.  17,  1962 

10. 

Alouette  (S-27) 

1004  KM 
540  N.MI. 
1.16  Re 

1029  KM 

554  N.MI. 
1.16  Re 

POLAR 

A 

105.5 

Sept.  29.  1962 

11. 

Explorer  XIV 

278  KM 

150  N.MI. 

99,000  KM 
53,000  N.MI. 

33° 

2185 

(36.4  hr.) 

Oct.  21,  1962 

12. 

1962  -  to* 
(STAR  AD) 

191  KM 

103  N.MI. 
1.03  R£ 

5550  KM 
3000  N.MI. 
1.87  Re 

71° 

Oct.  26,  1962 

13. 

ANNA  -  IB 

1090  KM 
584  N.MI. 

i.n  re 

1180  KM 

635  N.MI. 

119  Re 

50° 

Oct.  31,  1962 

14. 

Relay  -  I 

1321  KM 
713  N.MI. 
1.21  Re 

7439  KM 

4001  N.MI. 
2.17  Re 

47.5° 

185.1 

Dec.  13,  1962 

56 


TABLE  V-2 


Estimated  Time  in  Orbit  to  Degrade  Cells  to  75%  of  Initial  Short 

(Based  On  Telemetry  Data) 

Type  of  Mils  Shielding  Shield  Thickness 

Satellite  _ Cell _ and  Material _ gm/cm 


Circuit  Current 
T_ 


Satellite  #1 
(Ref.  26) 

10.  Aiouette  (S-27)  p/n 

11.  Explorer-XIV  p/n  Blue 


12  Glass 
Bare 


0.076 

0 


‘0.75 

Days 


Remarks 


1.  Explorer-Xl 
(Ref.  47) 

n/p 

Bare  0 

70 

Low  energy 
protons. 

2.  Midas  -  HI 
(Ref.  66) 

p/n  gridded 
p/n  gridded 
n/p 

40  Silica  0.275 

80  Silica  0.55 

60  Silica  0.41 

100 

250 

1000 

Data  corrected 
to  44°  C. 

3.  Explorer-XU 

p/n 

p/n 

3  Glass  0.02 

Bare  0 

1000 

1  Orbit 

Low  energy 
protons. 

4.  Midas  -  IV 
(Ref.  66) 

p/n 

p/n 

40  Silica  0.275 

80  Silica  0.55 

53 

110 

Data  at  87°  C. 

5.  TRAAC 
(Ref.  34) 

p/n 

6  Glass  0.038 

Knocked  out 
by  Starfish. 

6.  TRANSIT  -  IVB 
(Ref.  34) 

p/n 

6  Glass  0.038 

Knocked  out 
by  Starfish. 

7.  Ariel  -  1 

p/n 

• 

Knocked  out 
by  Starfish. 

8.  Telstar  -  I 

n/p 

30  Sapphire  0.3 

400 

(Partially) 

(Ref.  16) 

9.  Tetrahedral 
Research 

p/n  and  others 

Bare  0 

00 

Low  alt. 
orbit. 

30 

1  Orbit 


Low  energy 
protons. 


12.  1962  *  h  p/n 

(Ref.  33)  p/n 

n/p 

13.  ANNA  -  IB  p/n 

(Ref.  33)  p/n 

p/n 
n/p 

14.  Relay  *  1  p/n 

(Ref.  81  &  p/n 

82)  p/n 

n/p 
n/p 
n/p 


6  Glass 

0.038 

60  Glass 

0.38 

6  Glass 

0.038 

6  Glass 

0.038 

20  Silica 

0. 135 

30  Sapphire 

0.3 

30  Sapphire 

0.3 

Bare 

0 

30  SUica 

0.2 

60  SUica 

0.4 

Bare 

0 

30  SUica 

0.2 

60  SUica 

0.4 

4 

Affected  by 

100 

Russian  ex* 

50 

plosions. 

40 

Damage 

90 

occurs  in 

400 

South  Atlantic 

1400 

Anomaly 

0.04 

(Ref.  81) 

30 

(Ref.  81) 

46 

(Ref.  81) 

0.35 

(Ref.  81) 

160 

(Ref.  82) 

240 

(Ref.  82) 
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Reduction  In  Short  Circuit  Current  (Days) 


2 

0.2  gm/cm  .  Therefore,  the  data  from  Fig.  32  and  33  have  been  replotted  in  Fig. 
by  multiplying  the  p/n  life  values  by  the  factor  4.0.  It  may  be  seen  that  the  shape  of 
the  resulting  curve  of  the  predicted  life  of  n/p  silicon  solar  cells  vs.  shield  thickness 
appears  to  be  a  function  of  the  particular  satellite  orbit,  as  expected,  since  it  depends 
on  the  electron  and  proton  fluxes  and  energy  spectra  encountered  along  the  orbit.  The 
factor  of  4  difference  between  n/p  and  p/n  cells  is  only  expected  to  be  applicable  if 
most  of  the  damage  is  produced  by  protcns  above  about  2  Mev,  where  the  damage 
constants  are  found  w  differ  by  about  a  factor  of  four.  It  is  noted  that  the  factor  of 
four  difference  between  n/p  and  p/n  cells  does  not  apply  for  the  1962  BK  satellite 
with  thin  shielding  where  the  n/p  cells  appear  to  last  more  than  ten  times  as  long  as 
p/n  cells.  This  fact  may  be  due  to  one  or  both  of  two  factors:  (1)  The  p/n  cells 
degraded  rapidly  in  (our  days  at  a  time  when  the  flux  of  electrons  from  the  Russian 
explosions  was  high,  whereas  the  average  flux  during  the  50  day  life  of  the  n/p  cells 
may  have  been  lower,  and  (2)  even  if  the  average  flux  encountered  per  day  were 
constant  in  time,  but  the  damage  was  produced  mainly  by  electrons,  a  factor  of  the 
order  of  20  to  100  times  different  life  is  predicted  on  the  basis  of  the  relative  damage 
constants  for  1  Mev  electrons  and  a  factor  of  10  to  25  for  3  Mev  electrons.  Fig.  34 
may  be  used  as  a  flrat  approximation  to  estimate  the  life  cf  n/p  1  ohm -cm  silicon 
solar  cells  in  certain  orbits  which  have  orbital  parameters  fairly  close  to  those  for 
which  the  flight  test  data  were  obtained.  It  is  seen  that  the  most  severe  orbit  with 
Teapect  to  radiation  damage  is  that  of  Relay-1,  followed  in  order  by  1962  0  K,  Telstar 
a,  Midas -III  and  Anna-IB.  The  data  for  the  Teoahedral Research  Satellite  No.  1  is 
not  plotted  because  it  showed  no  measurable  degradation,  even  on  bare  p/n  cells, 
during  the  first  42  days,  which  indicates  that  the  damage  rate  to  solar  cells  is 
negligible  for  a  satellite  with  an  apogee  of  only  362  nautical  miles. 
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Tune  For  2S:'  Reduction  In  Short  Circuit  Current  (Days) 


VI  DESIGN  METHODS  FOR  SOLAR  CELL  POWER  SYSTEMS 


A.  General  Requirements 

Some  of  die  major  factors  which  need  to  be  considered  in  designing  a  solar  cell 
power  system  for  operation  in  a  radiation  environment  include: 

1.  The  severity  of  the  radiation  environment  expected  in  the  proposed  orbit. 

2.  The  availability  of  radiation  resistant  solar  cells  (e.g.  n/p  1  ohm -cm 
silicon  cells,  higher  resistivity  n/p  cells,  or  graded-base  cells). 

3.  The  thickness  of  shielding  which  can  be  used  within  the  weight  limit,  or  the 
increased  area  of  cells  which  can  be  used  to  produce  the  required  power 
within  the  area-limit.  (Generally  a  tradeoff  optimization  study  is  desirable.) 

4.  The  expected  range  of  operating  temperatures  for  the  solar  ceils,  which 
affects  the  selection  of  the  solar  cell  operating  point  (the  voltage  per  cell  or 
the  load  resistance). 

B.  Selection  of  Type  of  Solar  Cells 

Gent  'ally  if  a  proposed  satellite  is  to  spend  an  appreciable  fraction  of  its  time  in 
die  altitude  zone  between  about  1000  km  and  10,000  km,  it  will  be  necessary  to  use 
radiation  resistant  n/p  silicon  solar  cells  in  order  to  attain  reasonable  life .  In  other 
cases  there  may  be  an  economic  advantage  in  using  p/n  cells  unless  production  quan¬ 
tities  of  radiation  resistant  cells  become  available  at  equal  cost.  There  is  evidence  that 
n/p  silicon  cells  with  a  base  resistivity  greater  than  1  ohm-cm  and  also  graded -base 
(drift-field)  cells  are  more  radiation  resistant  than  1  ohm-cm  cells.  Therefore,  it  is 
expected  that  production  quantities  of  such  cells  will  become  available  for  use  in  extreme 
environments.  Solar  cells  of  gallium  arsenide,  cadmium  sulfide  and  other  materials  are 
still  in  the  experimental  research  stages. 

C.  Degradation  Rates  in  Particular  Orbits  a3  a  Function  of  Shield  Thickness 

For  certain  orbits  which  are  very  similar  to  those  of  the  satellites  discussed  in 
Section  (V,  an  approximation  of  the  time  for  silicon  cells  to  degrade  by  25%  in  short 
circuit  current  can  be  obtained  by  reference  to  Fig.  32,  33  and  34.  For  degradations 
greater  than  25%,  which  are  of  interest  in  overdesigning  the  total  cell  area,  it  may  be 
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assumed  that  the  short  circuit  current  degrades  by  approximately  20%  per  decade  in 
time .  This  is  a  good  approximation  only  if  the  radiation  dux  encountered  per  day  is 
constant,  or  if  the  time  for  25%  degradation  is  long  compared  to  the  period  for  apogee  pre¬ 
cession  (the  apogee  precession  period  for  Anna  -  IB  is  120  days).  If  one  desists  for  45% 
degradation  in  short  circuit  current,  the  life  will  be  extended  by  a  factor  of  10  as  com¬ 
pared  to  25%  degradation.  It  is  clear  that  the  use  of  area  overdesign  can  be  a  potent 
method  of  extending  solar  cell  life.  It  may  be  preferable  in  many  cases  to  increase  the 
number  of  solar  cells  used,  rather  than  increase  the  shield  thickness  in  order  to  obtain 
the  required  power  at  the  end  of  life . 

D.  Calculation  of  Solar  Cell  Degradation  Rates 

A  procedure  which  can  be  used  to  estimate  solar  cell  degradation  in  a  particular 
orbit  is  outlined  below.  The  procedure  given  for  calculating  electron  damage  is  only 
applicable  for  a  fission  electron  spectrum,  which  is  approximately  correct  for  orbits 
in  the  region  of  the  Mcllwain  L  parameter  less  than  2  (an  altitude  of  4000  miles  at  the 
magnetic  equator).  For  a  more  general  case  .changes  of  the  electron  energy  spectrum 
along  the  orbit  should  be  included  in  the  computation, 

a .  Select  a  value  of  the  initial  diffusion  length  of  the  solar  cells,  typically 
about  120  microns  (0.012  cm)  for  n/p  one  ohm -cm  silicon  cells. 

b .  Assume  an  allowable  degradation  for  the  short  circuit  current  or  maximum  power 
and  determine  the  final  value  of  diffusion  length  L  from  Fig.  15  or  19 

using  the  curve  corresponding  to  short  circuit  current  under  space  sunlight 
with  the  appropriate  value  of  L^. 

c  .  Calculate  the  allowable  change  in  (1/L^)  which  corresponds  to  the  final  short 

circuit  current: 

mi/l2)  .  rM02]  (i2> 

d  .  Calculation  of  Proton  Damage  Rate 

1.  Choose  an  energy  distribution  for  the  protons  which  is  applicable  to  the 
orbit  and  express  the  flux  per  unit  energy  interval  in  die  form: 

P(E>  =  Ioe-E/E° 
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This  equation  should  apply  over  the  range  from  E  ,  which  is  the  minimum 
proton  energy  capable  of  penetrating  the  assumed  shield  thickness  as  given 
by  Fig.  35  ,  up  to  approximately  E 2  =  200  Mev. 

2.  Use  a  proton  damage  coefficient  applicable  to  the  solar  cell  base  material 
as  shown  in  Fig.  11  or  12. 

3.  Calculate  the  average  damage  integral  per  unit  time: 

Dp  -  -f/  f2  Kp(E)„  (E.OdEd.  (11) 

0  Ex 

where  E,  is  the  minimum  energy  that  will  penetrate  the  shield  as  deter¬ 
mined  from  Fig.  35  .  The  integral  with  respect  to  time  may  be  approxi¬ 
mated  by  breaking  the  proton  belt  into  several  regions  in  B-L  coordinates 
and  calculating  the  time  spent  in  each  region,  taking  into  account  orbit 
regression  and  precession  of  the  orbital  apogee.  In  order  to  calculate 
the  average  damage  rate,  the  time  t  should  equal  one  period  of  preces¬ 
sion  of  the  apogee,  unless  it  is  desired  to  predict  the  degradation  more 
precisely  by  integrating  from  the  expected  initial  orbital  parameters. 

Calculation  of  Electron  Damage  Rate 

1.  Assume  an  energy  distribution  for  the  electrons  encountered  along  the 
orbit,  (e.g.  a  fission  spectrum  for  orbits  passing  through  the  artificial 
electron  belt  at  values  of  the  Mcllwain  L  parameter  less  than  2)  and 
calculate  the  average  omnidirectional  flux  of  fission  electrons  per  unit 
time  (epp). 

2.  Refer  to  Rosenzweig's  curve  of  equivalent  1  Mev  electron  flux  per  omni¬ 
directional  fission  electron  for  the  assumed  shield  thickness  and  solar 
cell  type  (Fig. 36  )  and  read  off  the  ratio  R  of  equivalent  1  Mev  integrated 
electron  flux  per  unit  omnidirectional  fission  electron  flux.  Calculate  the 
equivalent  average  1  Mev  flux  for  the  orbit: 


Read  the  K„ .  value  for  1  Mev  electrons  from  Fig.  8  for  the  type  of  solar 

C  1 

cell  base  material  used,  and  calculate  the  average  electron  damage  rate 


in  the  orbit: 


de  =  keivei 


This  has  the  units  of  change  in  1/L  per  unit  time, 
g .  Calculate  the  average  damage  rate: 


°T  =  DP  +  °E  (14) 

where  D  is  from  Step  D-3  and  D  is  from  Step  F. 

a  H 

h.  Solve  for  solar  cell  life  to  the  assumed  final  value  of  short  circuit  current 


from: 


A  (1/L  ) 


where  ^  (1/L2)  J  is  from  Step  C  and  D^,  from  Step  G. 

E.  Sample  Calculation  of  Radiation  Damage  in  Equatorial  Orbits 
1.  Calculation  of  Proton  Damage  - 

It  is  assumed  that  the  proton  spectral  flux  is  given  by  Mcllwain  and  Pizzella, 


(Ref.  100): 


where 


o(E)  =  IQe  E/E0 

Eq  =  306  L-5' 2  (Mev) 

L  =  Mcllwain’s  parameter 


For  a  satellite  at  an  altitude  of  2000  miles  (half  an  earth  radius)  L  =  1.5  and  Eq  =  37.3 

Mev.  The  value  of  I  is  found  by  normalizing  to  an  integral  flux  above  40  Mev.  of 
4  2” 

2  x  10  protons/cm  -  sec,  as  given  by  O'Brien  (Ref.  104): 

f  e-E/E0  dE  =  2  x  104  (18) 


or  Iq  =  1570  p/cm  -  sec  -  Mev 

The  radiation  damage  to  n/p  1  ohm -cm  silicon  solar  cells  covered  by  30  mils  of 
sapphire  shielding  on  the  front  and  with  infinite  back  shielding  will  be  estimated.  The 
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proton  damage  coefficient  is  assumed  to  be 

K  =  A/E 
P 

where  A  =  (1.5  x  10^)  for  E  <  60  Mev  (19) 

K  =  B  =  2.5xl0^  for  E  >  60  Mev 
P 

It  is  assumed  that  the  energy  spectrum  and  flux  of  protons  behind  the  shield  is  the  same 
as  on  the  outside  of  the  shield  except  that  the  low  energy  particles  are  stopped  by  the 
shield.  This  assumption  is  equivalent  to  assuming  that  the  lower  energy  and  therefore 
higher  damage  coefficient  for  those  protons  which  penetrate  the  shield  is  just  compen¬ 
sated  by  the  fact  that  an  omnidirectional  flux  of  incident  protons  is  strongly  absorbed 
at  small  angles  of  incidence  so  that  only  protons  within  a  certain  entrance  cone  can 
penetrate.  This  assumption  is  conservative  in  that  it  somewhat  over  estimates  the 
damage  produced  by  protons  which  have  energies  slightly  greater  than  that  required  to 
penetrate  the  shield.  (Compare  Rosenzweig's  calculated  results  as  shown  in  Fig.  37 
The  change  in  inverse  squared  diffusion  jengh  for  the  cells  is  calculated  from: 


00 

id/L2)  -  -y-  J  p(E)K(E)dE 


(20) 


where  T  =  time  of  exposure,  seconds 


Ec  =  cutoff  energy  =  14  Mev  for  30  mil  sapphire  (0.3  gm/cm  ) 


or 


*(1/L2)  =  J  AE  1e“fi/E0  dE  +J  Be^^O  dE 


“E/E, 


60 


(21) 


If  we  let  x  = 


Then 


where 


Ad/L.2)  -  Ar,6°  1  " 

*  X 

c 


f  6o  i  -x  r 

kJ  Te  dx+BEoJ 


e  dx  +  BE„  I  e  dx 
*60 


14 

Xc  =  E„ 


60 

*60  “  E„ 


(22) 


and 


J 


-x  .  -60/E. 

e  dx  =  e  0 


'60 


(23) 


t 
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The  first  integral  in  equation  22  may  be  evaluated  as  an  infinite  series: 


— x  2  3  4 

.  re  dx  .  x  x  //ii\ 

F(X)=J-J - ‘"’‘-x+rr2T  -3T37  +4T41  <«> 

A  plot  of  the  function  F(x)  is  shown  in  Fig.  38  .  Using  this  figure  we  evaluate  (with 

Eq  =  37.3): 

F(x,n)  ~F(xlA)  =  F(1.61) -F(0. 375)  =  -0.62  -(-1.30)  =  0.68  (25) 

OU  14 

Evaluating  the  damage  rate,  we  have: 

A(l/L2)  =  *0  |  1.02  x  10-5  +0.186  x  10-5!  =  6  x  lO”6^  (26) 

T  2  -  '  ° 

It  is  seen  that  the  damage  by  protons  above  60  Mev  is  only  18%  of  that  by  protons  from 

14  to  60  Mev  for  the  assumed  spectral  distribution. 

2 

For  Iq  =  1570  p/cm  -sec-Mev 

A(l/L^)  =  9.4  x  10  3  cm  2  sec  * 

T 

For  solar  cells  with  an  initial  diffusion  length  of  120u,  a  decrease  of  25%  short  circuit 
current  is  obtained  when  L  is  about  15. 5u.  The  time  required  to  degrade  to  L  =  15.5u 
by  protons  only  is: 

— l—  ~  108 

T  =  (15.5)  (120) 2 

P  9.4  x  10~3 

•p  -  4. 10  x  IQ3  =  4.35  x  10  7  sec  =  505  days 
P  9.4  xlO-3 

This  is  the  calculated  life  under  proton  damage  to  reach  75%  of  initial  current  for  n/p 
1  ohm -cm  cells  with  30  mils  of  sapphire  shielding  in  a  circular  equatorial  orbit  at  2000 
miles  altitude. 

2.  Calculation  of  Electron  Damage  - 

For  equatorial  orbits  between  L  =  1.3  and  1.8  (1200  to  3200  miles  altitude) 

it  is  assumed  that  a  fission  electron  spectrum  applies  with  an  omnidirectional  flux  of 
8  2 

5  x  10  e/cm  sec,  which  is  one  half  the  value  given  by  Hess  for  a  time  of  one  week  after 
Starfish.  (This  is  probably  too  high  a  flux  and  too  hard  a  spectrum  at  L=1.8, 
based  on  personal  communication  from  Walter  L .  Brown .) 


From  Rosenzweig's  data  (Fig.36  )  the  conversion  factor  from  omnidirec¬ 
tional  fission  electrons  to  equivalent  normal  incident  1  Mev  electrons  is  0.37  for  a 

2 

shield  thickness  of  0.3  gm/cm  on  n/p  cells. 

2 

The  change  in  1/L  per  1  Mev  electron  is  found  to  be  approximately 


4.1 x  10“ 
2.5  x  101 


=  1.64  x  10 


(Fig. 8  shows  K_.,  =  1.7  x  10  at  1  Mev) 
c  1 

The  average  electron  damage  rate  is: 


DE  *  KE1  »El 


—10  ft 

(1.64  x  10  )  (0.37)  (5  x  10  ) 


—  -2  —1 
D„  =  0.0303  cm  sec 
E 

The  time  to  degrade  the  cells  by  25%,  considering  only  electrons  is: 


(i/L2)  =  4.1  x 


0.0303 


=  1 . 35  x  10  sec  =  156  days 


3.  Calculation  of  Combined  Effect  of  Electrons  and  Protons  - 


The  total  damage  rate  is 


D  =  D  +  D_ 
T  E  P 


from  which 


6(1/1/ 

T_ 


J_  +  _L 

T  T 
E  lP 


T  VjL 

T  "T2  +  TP 


Ll/L2)  ,  M1/L2) 


The  life  under  combined  electrons  and  protons  is  then 


t  =  =  il9  davs 

T  (156  +505)  y 

This  is  the  calculated  life  to  reach  75%  of  initial  current  for  n/p  1  ohm -cm 
silicon  cells  with  30  mils  of  sapphire  shielding  in  a  circular  equatorial  orbit  at  2000 
miles  altitude. 
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Similar  calculations  have  been  made  for  other  altitudes  using  the  proton  flux 
map  shown  in  Fig.  30  and  the  results  are  plotted  in  Fig.  39.  It  is  seen  that  the 
damage  rate  by  protons  increases  with  altitude  and  has  not  reached  its  maximum  at 
an  altitude  of  5000  km,  above  which  the  absolute  value  of  proton  fluxes  are  not  well 
defined.  This  effect  is  due  to  the  increasing  fluxes  of  protons  with  energies  slightly 
above  the  14  Mev  shield  cutoff  energy  which  are  encountered  with  increasing  values  of 
the  Mcllwain  L  parameter. 

The  solar  cell  life  of  80  days  calculated  for  an  altitude  of  5100  km  may  be  com¬ 
pared  with  the  observed  life  of  190  days  interpolated  from  Relay- 1  flight  test  data. 

One  may  conclude  that  the  Relay- 1  damage  rate  is  equivalent  to  spending  about  40%  of 
its  time  at  the  most  damaging  point  calculated,  namely  at  5100  km  altitude  at  the 
equator.  Actually  Relay  spends  considerable  time  at  larger  L  values  than  1.8  where 
it  encounters  many  low  energy  protons  and  it  also  passes  through  the  artificial  electron 
belt.  Therefore,  a  computer  calculation  would  be  required  after  the  fluxes  and  spectra 
are  better  defined  in  order  to  compare  predicted  with  observed  degradation  rates. 

F.  Comparison  of  Calculated  and  Observed  Damage  on  S-27  Alouette 

The  Alouette  is  in  approximately  a  1000  km  circular  polar  orbit  for  which  Hess 

12 

(Ref.  39)  has  calculated  an  omnidirectional  flux  of  fission  electrons  of  2  x  10  per 

day  at  a  time  of  one  week  after  Starfish.  Alouette  was  actually  launched  on  Sept.  29, 

1962,  nearly  three  months  after  Starfish.  For  blue-shifted  p/n  solar  cells  with 
2 

12  mils  (0.076  gm/cm  )  of  glass  shielding,  Rosenzweig’s  conversion  ratio  from 

12  2 

Fig.  36  is  found  to  be  0.37,  so  the  equivalent  1  Mev  flux  is  0.74  x  10  e/cm  day. 

The  damage  coefficient  for  p/n  cells  at  1  Mev  is  selected  from  Fig.  9  to  be 
_8  2 

1.22  x  10  and  the  change  in  1/L  to  produce  25%  reduction  in  short  circuit  current 
5  -2 

is  taken  as  4. 1  x  10  cm  .  Therefore,  the  time  to  reach  25%  reduction  in  short 
circuit  current  under  electron  damage  is 


.  aq/L2) 

E  Kcp 


T 


E 


4. 1  x  105 _ 

(1.22  x  10-8)  (0.74  x  1012) 


45  days 


This  is  probably  an  underestimate  of  "electron  life"  because  of  decay  of  the  artificial 
belt. 


0 


0 


0 


0 


0 


0 


74 


Time  Vox  \2S% 


t 

Flight  test  data  from  Alouette  showed  a  25%  degradation  after  about  30  days. 

The  added  damage  sustained  may  be  accounted  for  by  protons  encountered  in  the  South 

Atlantic  Anomaly.  _ 
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APPENDIX  A 

RADIATION  EFFECTS  ON  SOLAR  CELL  COVER  SLIDE  MATERIALS  AND  ADHESIVES 


A .  Introduction 

The  prediction  of  solar  cell  power  system  life  should  take  into  account  the 
degradation  of  the  optical  transmissive  properties  of  the  cover  materials .  The  covering 
materials  are  generally  a  transparent  cover  slide  and  a  transparent  adhesive .  The 
cover  slide  is  utilized  to  shield  against  electrons  and  protons,  and  to  carry  a  spectrally 
selective  filter .  It  also  protects  against  contamination  from  being  handled  on  the 
ground .  Investigations  have  been  conducted  with  protons  and  electrons  on  various  cover 
slide  materials  and  adhesives .  From  these  studies  an  indication  of  the  degradation 
that  can  be  expected  in  service  may  be  derived . 

B.  Damage  Mechanism 

Cover  materials  are  susceptible  to  "F"  center  formation  from  high  energy 
particle  fluxes .  The  formation  of  these  color  absorbing  centers  is  a  function  of  the 
impurities  in  the  material .  Most  color  centers  can  be  annealed  out  of  the  material 
by  heat ,  infrared  radiation  or  visible  light . 

In  general ,  the  sensitivity  to  coloration  is  not  fixed  for  a  given  type  of  engineering 
material  but  varies  according  to  the  impurities ,  lattice  imperfections ,  strain ,  stoichiometry 
and  the  prior  deformation  of  the  material  (Ref.  118).  When  irradiation  is  at  elevated 
temperatures  an  equilibrium  density  of  color  centers  is  attained  from  the  production  and 
destruction  of  the  color  centers  . 

Relatively  little  work  has  been  done  on  coloration  or  the  coloration  equilibrium 
from  simultaneous  energetic  protons  ,  electrons  and  ultraviolet  radiation  below  3000 
Angstroms  in  the  presence  of  the  longer  wavelength  components  of  sunlight.  However, 
the  change  in  transmittance  of  various  materials  from  proton  and  electron  radiation 
has  been  measured  by  F  .  A .  Campbell  of  the  Naval  Research  Laboratory  (Ref.  H2) 
and  by  J  ,  L  .  Patterson  and  G .  A .  Haynes  (Ref.  120)  of  the  Langley  Research  Center . 

These  data  give  an  indication  of  the  amount  of  damage  that  can  be  expected  in  service . 


C  .  Cover  Slide  Requirements 


The  most  desirable  materials  for  solar  cell  cover  slides  are  those  which  show 
the  least  amount  of  transmission  degradation  within  the  solar  cell  spectral  response 
band  upon  exposure  to  electron  and  proton  radiation .  This  spectral  band  is  typically 
0.35  to  1 .1  microns,  with  a  peak  near  0.8  microns  .  When  the  spectral  distribution 
of  space  sunlight  is  taken  into  consideration,  the  peak  energy  point  occurs  at  approxi¬ 
mately  0 .6  microns . 

D.  Electron  and  Proton  Damage  Experiments 

Assessment  of  electron  damage  has  been  made  by  F.  A.  Campbell  (Ref.  112 

and  113)  of  the  Naval  Research  Laboratory.  Electrons  of  1  Mev  energy  at  a  total  dose 
16  2 

of  10  e/cm  were  utilized .  Spectral  transmittance  v/as  measured  with  a  Beckman 

IR-4U  double  beam  spectrophotometer.  The  electron  exposure  was  conducted  in  air, 

directly  under  the  beam  of  a  Van  de  Graaff .  Transmittance  values  for  various  materials 

between  0  .4  and  0 .7  microns  are  summarized  in  Table  A-l  (Ref .  113) . 

Patterson  and  Haynes  (Ref.  120)  utilized  a  Radiation  Dynamics  No.  P.E  .A. -10 

2 

accelerator  with  a  constant  beam  current  density  of  0.03  microamp/cm  .  Electron 

energy  of  1 .2  Mev  was  used  and  most  of  the  samples  were  exposed  to  a  total  dose  of 
15  2  c 

2  .7  x  10  e/cm  .  A  2900  K  filament  temperature  tungsten  lamp  was  used  as  the 
source  and  the  transmission  efficiency  was  recorded  as  a  percent  change  in  broad 
spectral  band  light  passing  through  the  sample  to  a  typical  solar  cell .  The  measure¬ 
ment  accuracy  was  1/2%. 

Measurements  of  changes  in  spectral  transmission  after  proton  irradiation 

were  conducted  it  the  Naval  Research  Laboratory.  The  irradiation  consisted  of  a 

’  •  2 

single  dose  of  4  x  10  p/cr.-.*-  at  an  energy  of  4  .6  Mev .  The  samples  were  irradiated 

-4 

in  vacuum  at  a  pressure  of  10  mm  Hg  and  the  results  are  summarized  in  Table  A-2. 
The  test  data  show  a  superiority  for  the  Corning  7940  specimens  as  compared  to 
Coming  0211 . 

1 .  Silica 

Two  major  classes  of  silica  are  synthetic  qu^.-tz  (or  fused  silica)  and  vapor 
deposited  silica  .  Both  electron  and  proton  damage  is  reported  for  variour  types  and 
manufacturers  of  fused  silica  as  shown  in  Tables  A-l,  A-2  and  A-3.  Electron  radiation 


A2  •• 


I 


88.5  88.5  88.0  89.0 


TABLE  A-l  cont’d 
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TABLE  A-3 

Effects  of  1 .2  Mev  Electron  Radiation  on  Transparent  Materials  (From  Patterson  &  Haynes ,  Ref .  120) 
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damage  curves  for  Coming  0211  (microsheet)  with  various  types  of  filters  and  anti¬ 
reflection  coatings  are  plotted  separately  in  Fig.  A-l,  A-2,  and  A-3  between  wave¬ 
lengths  of  0 .35  to  1 .20  microns  .  The  spectral  degradation  of  transmittance  for 
electron  and  proton  irradiation  of  Coming  0211  is  predominately  at  the  shorter  wave¬ 
lengths  and  decreases  gradually  with  increasing  wavelength .  The  "blue-red"  filter 
specimen  however  is  an  exception  and  in  Figure  A-3  shows  appreciable  degradation 
at  all  wavelengths  up  to  about  1 .0  micron .  The  microsheet  (Coming  0211)  seems  to 

approach  a  saturation  limit  since  it  was  found  that  the  transmittance  at  all  wavelengths 

i  6  2  2 

was  only  slightly  lower  at  10  e/cm  than  at  the  10  e/cm  level . 

Electron  radiation  damage  for  Coming  7940  fused  silica  between  0.40  and  0.70 
microns  is  summarized  in  Table  A-l ,  numbers  5  and  6 .  The  7940  specimen  showed 
2.2%  or  less  change  in  transmittance  for  both  types  of  particle  irradiation  and  was  the 
best  type  of  cover  material  reported  . 

Table  A-3  summarizes  the  decrease  in  transmittance  oi  fused  silica  with  1 .2  Mev 

15  2 

electrons  at  a  total  dose  of  2  .7  x  10  e/cm  .  Coming  7940,  Suprasil  II,  Spectrasil  A , 
and  Dynasil  Optical  Grade  had  no  detectable  degradation  and  in  each  case,  were  made  by 
a  vapor  deposition  process  .  An  example  of  the  scatter  in  transmittance  values  as  a 
function  of  purity  of  silica  as  shown  in  Fig .  A-4  .  The  poorest  silica  had  50  to  70  parts 
per  million  of  impurities  and  showed  close  to  a  40%  decrease  in  transmission  .  The 
vapor-deposited  silica  had  less  than  1  part  per  million  of  impurity  and  showed  a 
superiority  over  the  fused  crystalline  silica .  Those  silica  samples  that  darkened  in 
Table  A-3  were  made  by  the  process  of  fusing  crystalline  quartz . 

In  general ,  in  each  case  in  Table  A-3  where  darkening  was  noted,  the  curves 
indicated  a  broad  absorption  band  centered  at  approximately  5500  A°.  Therefore,  if 
space  sunlight  illumination  had  been  used  with  a  solar  cell,  the  decrease  of  useful  light 
transmission  would  have  been  greater  than  indicated. 

Commercial-grade  fused  silica,  1/16  inch  thick,  from  the  Amersil  Division  of 
Engelhard  Industries  ,  Inc  .  was  irradiated  with  protons  in  conjunction  with  a  solar  cell 

experiment  at  Langley  Research  Center  (Ref.  116).  There  was  no  darkening  of  the  silica 

11  2  13  2 

with  a  5  x  10  p/cm  dose  of  240  mev  protons ,  but  with  a  2  x  10  p/cm  dose  of  22  Mev 

protons  the  cover  slide  darkening  reduced  the  short  circuit  current  by  3.7%  for  one  sample 
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and  7 .6%  for  a  second .  The  solar  cells  were  8%  efficiency  p/n  cells  and  the  light 
source  was  a  lamp  with  a  tungsten  filament  operating  at  2900°K  with  a  1- inch -thick 
water  filter . 

Tests  of  ordinary  quartz  shielding  by  Kallander  and  Weller  (Ref.  119)  with 
15  2 

an  integrated  flux  of  6  x  10  1  Mev  electrons  per  cm  also  showed  considerable 

darkening,  which  indicates  that  impure  silica  is  to  be  avoided . 

2.  GlasB 

Typical  varieties  of  glass  are  very  susceptible  to  radiation  damage  as  shown 

in  Table  A-3 .  The  ordinary  soda  lime  plate  glass  degraded  in  transmittance  by  26 .0% 

15  2 

with  a  dose  of  2 .7  x  10  1 .2  Mev  electrons/cm  .  Some  protection  is  available  with 

the  addition  of  trivalent  cerium  which  tends  to  neutralize  the  color  centers ,  thereby 
reducing  the  absorption . 

3.  Sapphire 

Artificial  sapphire  has  been  used  as  a  cover  shield  material  for  solar  cells 

(e  .g .  on  Telstars  1  and  II) .  The  Langley  Research  Center  irradiated  a  cover  shield 

made  of  artificial  sapphire  from  the  Linde  Company ,  with  1 .2  Mev  electrons  at  a  total 
15  2 

dose  of  2 .7  x  10  e/cm  .  With  this  flux ,  no  decrease  in  broad  band  transmittance  was 
observed .  This  material  as  shown  in  Table  A-3  was  the  only  material  other  than 
vapor-deposited  silica  that  suffered  no  measurable  decrease  in  transmittance . 

4.  Adhesives 

Adhesives  are  often  utilized  to  bond  the  cover  shield  to  the  face  of  the  solar 

cell .  The  adhesive  requires  dimensional  stability ,  transparency ,  high  electrical 

resistivity ,  long  term  adhesion  and  chemical  compatibility  with  the  solar  cell .  The 

primary  concern  to  date  has  been  'vith  the  change  in  transparency  under  proton,  electron 

and  ultraviolet  radiation .  The  effect  of  Cobalt  60  gamma  radiation  on  the  resistivity 

of  an  Epoxy  adhesive  is  reported  by  Kallander  and  Weller  (Ref.  119).  The  surface 

9  8 

resistivity  decreased  to  10  ohm-cm  at  a  dose  of  2  x  10  Roentgens . 

F  .  A .  Campbell  of  the  Naval  Research  Laboratory  has  measured  the  change 
in  spectral  transmission  of  adhesives  Dow  DER-332,  Epocast  15-E  and  Spectrolab  ES-10 
with  electron  irradiation .  The  adhesives  were  cast  1-2  mils  thick  in  the  same  method 
as  used  to  cement  cover  slides  to  solar  cells  .  The  adhesives  ES-10  and  15-E  were  cast 
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under  a  6  mil  Corning  7940  fused  silica  cover  and  the  DER-332  under  a  66  mil  fused 
silica  cover .  The  additional  shielding  did  not  help  the  DER-332  and  it  degraded  more 
than  the  ES-10  adhesive  which  was  less  well  protected .  The  data  are  tabulated  in 
Table  A-l,  numb.  10,11  and  12.  The  ES-10  shows  a  marked  superiority  and  did  not 
degrade  more  than  2%.  The  1S-E  showed  a  decrease  in  transmission  with  each  successive 
electron  exposure ,  particularly  at  the  shorter  wavelengths .  Crazing  and  delamination 
were  also  visible  in  this  specimen  which  would  diffuse  the  spectrometer  light,  there¬ 
fore  giving  a  lower  transmissivity  reading.  However,  this  diffused  light  may  still 
be  utilized  by  a  solar  cell . 

E  .  Conclusion 

The  materials  tested  which  indicated  a  superiority  by  showing  minimum 
darkening  with  energetic  electrons  are  Corning  7940  of  Corning  Glass  Works , 

Suprasil  II  and  Optical  Grade  Silica  of  Engelhard  Industries,  Inc . ,  Spectrasil  A 
of  Thermal  American  Fused  Quartz  Co . ,  GE  104  of  General  Electric  Co . ,  Dynasil 
(Optical  Grade)  of  Dynasil  Corporation ,  and  synthetic  sapphire  of  the  Linde  Company , 

The  material  tested  with  the  least  change  in  transmittance  for  4 .6  Mev  protons  was 
Coming  7940 .  The  adhesive  tested  with  the  least  change  in  transmittance  under 
electron  irradiation  was  the  Spectrolab  ES-10  specimen. 


APPENDIX  B 

SPACE  RADIATION  EFFECTS  ON  TRANSISTORS  8t  DIODES 

A .  Introduction 

In  addition  to  solar  cells ,  the  semiconductor  devices  of  primary  concern  with 
respect  to  radiation  damage  are  transistors  and  diodes ,  which  also  depend  on  diffusion 
of  minority  carriers  for  their  operation .  Primary  consideration  will  be  given  to  those 
types  of  transistors  and  diodes  associated  with  solar  cell  power  systems .  In  some 
applications ,  power  regulation  components ,  such  as  diodes ,  Zener  diodes  and  transistors , 
have  been  located  on  the  solar  panels  or  satellite  skin  where  they  are  exposed  to  space 
radiation  with  a  minimum  of  shielding.  Even  when  located  inside  a  satellite ,  it  is  found 
that  transistors  and  diodes  having  gas -filled  cases ,  can  be  seriously  affected . 

B.  Radiation  Damage  Mechanisms 

The  primary  radiation  effects  to  transistors  and  diodes  which  are  of  concern  in 
satellites  are  ( 1)  the  surface  effects  associated  with  ionization  of  gas  in  the  device 
enclosure,  and  (2)  the  permanent  damage  to  the  minority  carrier  diffusion  length 
produced  by  introduction  of  lattice  defects .  The  effect  of  injection  of  electron-hole 
pairs  by  ionization  in  the  semiconductor  material  is  of  major  significance  only  tinder 
extremely  high  flux  conditions ,  such  as  the  transient  neutron  and  gamma  radiation 
produced  by  a  nuclear  explosion ,  and  therefore  will  not  be  discussed  here . 

C.  Surface  Effects  of  Radiation 

For  transistors  and  diodes  which  have  a  gas -filled  enclosure,  the  rather  complex 
effects  of  ionization  cf  the  gas  interacting  with  the  surface  of  the  semiconductor  can 
produce  changes  i  1  performance  parameters  at  radiation  doses  as  low  as  10*  Roentgens . 
Such  a  dose  may  be  accumulated  in  a  period  of  hours  to  months  in  the  magnetosphere , 
depending  on  the  amount  of  shielding .  Therefore  it  is  important  to  select  transistors 
which  are  least  sensitive  to  this  effect . 

Reference  should  be  made  to  the  papers  by  Peck ,  Blair ,  Brown  and  Smits  (Ref .  141 
or  142)  for  a  detailed  discussion  of  this  surface  effect  which  occurs  in  transistors  under 
a  combination  of  a  radiation  environment,  a  bias  voltage,  and  a  gas  in  contact  with  the 


surface  .  The  effects  observed  depend  on  the  junction  bias  ,  the  envelope  voltage,  and  in 

many  respects  on  the  total  radiation  dose  ,  rather  than  dose  rate .  Upon  removal  of 

the  bias  voltage  or  the  radiation  field ,  a  device  tends  to  recover . 

It  has  been  found  that  large  variations  in  test  results  are  observed  between 

units  of  the  same  type  made  by  the  same  manufacturer . 

It  was  found  that  some  devices  can  be  screened  for  sensitivity  to  radiation  by 

means  of  a  short-time ,  high  level  radiation  dose ,  with  correlation  to  the  subsequent 

behavior  under  low-level  radiation  exposure. 

The  data  generally  indicate  that  transistors  and  diodes  in  evacuated  cases  and 

planar  units  (with  protection  over  the  semiconductor  surface)  are  least  affected  by 

ionizing  radiation  (Ref.  143).  By  the  use  of  evacuated  cases,  the  dose  at  which  a 

noticeable  increase  in  the  collector  reverse  current  (I  )  occurs  can  be  extended  to 
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the  range  of  10  to  10  Roentgens  as  seen  in  Fig.  B-l  (from  Ref.  142). 

D  .  Degradation  of  Diffusion  Length  in  Transistors 

Of  the  changes  in  semiconductor  properties  affected  by  radiation,  changes 
in  minority  carrier  diffusion  length  (or  lifetime)  are  the  predominant  cause  of  permanent 
damage  in  transistors  and  diodes  .  A  large  amount  of  testing  has  been  conducting  on 
germanium  and  silicon  devices  in  neutron  and  gamma  ray  fluxes  and  damage  coefficients 
have  been  determined.  However,  there  is  no  clear  correlation  which  permits  the  use 
of  data  determined  in  neutron  or  gamma  ray  fluxes  to  predict  damage  by  protons  or 
electrons  (Ref .  125). 

The  transistor  parameter  which  is  affected  most  seriously  by  a  degradation  of 

the  diffusion  length  is  the  common  emitter  forward  gain(h_D),  which  is  also  called  B. 

FE 

Bilinski  et  al  (Ref.  125)  give  the  following  equations  for  correlating  transistor  performance: 
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where  h__„  .  , 

FEO  =  initial  current  gain 

hp^  =  current  gain  after  bombardment 

t  =  overage  transit  time  of  minority  carriers  in  base 
0  =  integrated  flux 

=  lifetime  damage  constant 
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Integrated  Gamma  Dose  (R) 


(Radiation  Dose  Rate  8 .5  x  10^  rads/hr) 

FIG  .  B-l  EFFECT  OF  GAMMA  RADIATION  ON  COLLECTOR  REVERSE 
CURRENT  OF  EVACUATED  AND  GAS-FILLED  SILICON 
TRANSISTORS 
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and  K  is  given  by 
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The  average  transit  time  of  minority  carriers  in  the  base  region  of  the  transistor 
is  given  by: 


For  a  uniform  base 
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For  a  linear  graded  base 
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where  W  =  base  width 

L  =  minority  carrier  diffusion  Length 
f^k  =  alpha  cutoff  frequency 

The  above  expressions  agree  well  with  experiment ,  except  in  the  case  of  P-N-P 
germanium  transistors  where  the  effect  of  depletion  layer  widening  of  the  collector 
region  is  appreciable  .  (This  case  has  been  treated  by  Easley  and  Dooley,  Ref.  129). 
Except  in  this  case,  the  lifetime  damage  constant  for  germanium  and  silicon  transistors 


is  given  by: 
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K  is  related  to  the  diffusion  length  damage  constant  K.  by  K  -  DK 

T  T 


where  D  is  the  diffusion  coefficient  for  minority  carriers . 

Measurements  of  germanium  and  silicon  transistor  performance  after  proton 
irradiation  at  22  and  240  Mev  have  been  reported  by  Hulten,  Honaker  and  Patterson  of 
the  NASA  Langley  Research  Center  (Ref .  42).  Tests  with  40  aad  440  Mev  protons  are 
reported  by  Honaker  and  Bryant  (Ref.  132). 

An  attempt  to  correlate  the  Langley  data  on  transistors  in  terms  of  lifetime 
damage  coefficients  for  the  base  materials  was  made  by  Bilinski,  et  al(Ref.  125). 
However  ,  the  data  showed  considerable  scatter ,  some  values  for  p-type  silicon  being 
an  order  of  magnitude  smaller  than  the  damage  coefficients  measured  for  p-type  silicon 
in  solar  cells  by  Bell  Telephone  Laboratories  . 
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It  is  apparent  that  further  careful  research  is  required  to  obtain  quantitative 
data  on  which  to  base  predictions  of  the  degradation  of  transistor  performance  under 
penetrating  proton  radiation  .  However ,  general  rules  and  guidelines  for  the  selection 
of  transistors  can  be  derived  from  available  test  data  . 

E  .  Selection  of  Transistor  Type 

There  is  evidence  that  germanium  transistors  are  more  resistant  to  proton 
(and  neutron)  radiation  than  are  silicon  transistors  which  are  designed  for  the  same 
alpha  cutoff  frequency.  Also  the  sensitivity  to  radiation  damage  is  much  less  for 
transistors  with  a  narrow  base  width  and  a  high  alpha  cutoff  frequency .  There  is  some 
evidence  that  germanium  P-N-P  types  are  more  resistant  than  N-P-N  types  .  (This  is 
found  to  be  true  under  neutron  radiation  ,  Ref.  124) .  Therefore  one  may  conclude  that 
the  most  radiation  resistant  transistor  should  be  a  P-N-P  germanium  type  with  a  thin 
base  region  and  correspondingly  high  alpna  cutoff  frequency .  In  order  to  minimize 
surface  effects  ,  the  transistor  should  be  vacuum  encapsulated.  If  the  temperature 
environment  dictates  the  use  of  silicon ,  the  cutoff  frequency  should  be  high  .  There 
does  not  appear  to  be  a  clear  distinction  between  the  resistance  of  P-N-P  and  N-P-N 
silicon  types  ,  although ,  on  the  basis  of  solar  cell  experience  ,  one  would  expect  the 
N-P-N  type  to  be  superior  under  proton  and  electron  irradiation  . 

Unless  care  is  taken  to  avoid  surface  effects ,  a  transistor  may  fail  by  increase 
in  collector  reverse  current  at  a  total  radiation  exposure  much  lower  than  that  required 
to  degrade  the  current  gain  . 

F .  Diodes 

The  behavior  of  diodes  and  rectifiers  ,  and  Zener  reference  diodes  in  neutron 
and  gamma  ray  fields  has  been  reviewed  by  Reid  (Ref.  64).  In  general,  there  was 
often  an  initial  change  in  'he  reverse  characteristic  of  diodes  (which  was  probably 
attributable  to  surface  effects) ,  followed  by  a  degradation  of  the  forward  characteristic  . 

The  life  determined  by  die  allowable  increase  in  forward  impedance  can  be  greatly 
increased  by  using  a  diode  of  narrow  base  width  .  A  factor  of  10  reduction  in  base  width 
yields  an  increase  in  radiation  tolerance  by  a  factor  of  about  100. 

In  comparing  alloy  and  diffused  diode  structures  ,  Huth  (Ref .  133)  observes  that 
the  diffused  structure  marked’y  improves  radiation  behavior.  Reid  (Ref.  64)  concludes 
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that  the  optimized  configuration  of  a  silicon  diode  for  use  in  radiation  fields  would 
consist  of  a  diffused  p  -  i  -  n  structure  with  the  minimum  base  thickness  consistent 
with  peak  inverse  -  voltage  requirements  .  Based  on  present  knowledge  of  surface 
effects  ,  it  appears  that  vacuum  encapsulation  would  be  desirable  . 

G  .  Shielding  Against  Space  Radiation 

In  estimating  the  life  of  transistors  and  diodes  in  a  space  environment,  it  must 
be  recognized  that  surface  effects  are  dependent  on  the  ionization  rate  produced  by 
protons  ,  electrons  and  Bremstrahlung  (x-rays)  in  a  gas  filled  device  enclosure ,  whereas 
the  degradation  of  minority  carrier  lifetime  depends  on  the  flux  and  energy  spectrum  of 
those  penetrating  electrons  and  protons  which  have  sufficient  energy  to  produce  lattice 
displacements  (>  145  Kev  for  electrons  and  >  98  ev  for  protons  in  silicon) . 

The  severity  of  the  surface  effect  problem  is  critically  dependent  on  the  flux 
and  energy  spectrum  of  the  electrons  in  the  artificial  belt  which  are  predominant  over 
protons  in  determining  the  ionization  rate  (RADS/day)  in  gas  filled  devices  (Ref.  143). 
Fortunately,  there  is  evidence  that  the  electron  energy  spectrum  is  now  much  softer 
(relatively  fewer  high  energy  electrons)  than  a  fission  spectrum ,  at  least  tn  the  upper 
regions  of  the  artificial  belt  (Ref  .  126) .  However ,  it  appears  desirable  to  use  evacuated 
devices  where  possible  and  to  locate  any  gas-filled  devices  inside  a  satellite  to  take 
maximum  advantage  of  available  shielding,  if  operation  is  required  in  the  artificial 
electron  belt .  Gas-filled  devices  should  be  carefully  selected  and  tested  for  sensitivity 
to  surface  effects  . 

If  transistors  or  diodes  are  used  outside  the  satellite  shell,  they  should  be 

evacuated  and  have  protection  by  at  least  the  same  thickness  of  shielding  as  the  solar 

2 

cells  For  example  ,  at  least  0.3  grams/cm  of  shielding (45  mils  of  aluminum)  is 
desirable  for  a  satellite  orbit  like  Telstar  or  Relay  which  passes  through  the  inner 
proton  and  electron  belts  . 
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